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ABSTRACT 

Performing two-dimensional hydrodynamic simulations including a detailed treatment of the equation of state of the stellar plasma 
and for the neutrino transport and interactions, we investigate here the interplay between different kinds of non-radial hydrodynamic 
instabilities that can play a role during the postbounce accretion phase of collapsing stellar cores. The convective mode of instability, 
which is driven by the negative entropy gradients caused by neutrino heating or by variations in the shock strength in transient phases 
of shock expansion and contraction, can be identified clearly by the development of typical Rayleigh-Taylor mushrooms. However, in 
those cases where the gas in the postshock region is rapidly advected towards the gain radius, the growth of such a buoyancy instability 
can be suppressed. In this situation the shock and postshock flow can nevertheless develop non-radial asymmetry with an oscillatory 
growth in the amplitude. This phenomenon has been termed "standing (or spherical) accretion shock instability" (SASI). It is shown 
here that the SASI oscillations can trigger convective instability, and like the latter, they lead to an increase in the average shock radius 
and in the mass of the gain layer. Both hydrodynamic instabilities in combination stretch the advection time of matter accreted through 
the neutrino-heating layer and thus enhance the neutrino energy deposition in support of the neutrino-driven explosion mechanism. A 
rapidly contracting and more compact nascent neutron star turns out to be favorable for explosions, because the accretion luminosity 
and neutrino heating are greater and the growth rate of the SASI is higher. Moreover, we show that the oscillation period of the SASI 
observed in our simulations agrees with the one estimated for the advective-acoustic cycle (AAC), in which perturbations are carried 
by the accretion flow from the shock to the neutron star and pressure waves close an amplifying global feedback loop. A variety of 
other features in our models, as well as differences in their behavior, can also be understood on the basis of the AAC hypothesis. The 
interpretation of the SASI in our simulations as a purely acoustic phenomenon, however, appears difficult. 
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1. Introduction 

Hydrodynamic instabilities play an important role in core- 
collapse supernovae, because on the one hand they may be 
crucial for starting the explosion and on the other hand 
they may provide a possible explanation for the observed 
anisotropy of supernovae. There is a growing consensus that 
the neutrino-driven explosion mechanism of core-collapse su- 
pernovae does not work in spherical symmetry for progen- 
itors more massive than about 10 M Q . None of the recent 
simulations with one-dimensional (ID) hydrodynamics and a 
state-of-the-art description of the neutrino transport develops 
an explosion (Rampp & Janka 2002; Liebendorferet al. 2001; 
Liebendorfer etaLl 120051; iThompson et al.1 120031; iBuras et all 
20031 l2006allbn T However, multi-dimensional effects were recog- 
nised to be helpful. In particular it was shown that convec- 
tion is able to develop below the stalled supernova shock and 
that it can increase the efficiency of neutrino heating signifi- 
cantly dHerant et al.lll994t iBurrows et alii 19951; Uanka & Mtuleil 
1995l[l996l) . Current two-dimensional (2D) simulations are thus 
considerably closer to the explosion threshold than ID models 
(Buras et al. 2003, 2006a, b), and shock revival and the onset of 



an explosion has been re ported recently for a 2D calculation with 
an 11.2M Q progenitor dBuras et"aTll200 6b). In earlier 2D sim- 
ulations, in which the angular size of the numerical grid was 
constrained to less than 180° and in simulations in which the 
approximative description of the neutrino transport resulted in a 
fast onset of the explosion, convection was domi nated by rather 
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small angu lar scales of several ten degrees (Janka & Midler 
1994, 1996). However, i n recent 2D calculations oflBuras et al. 
(200 6allbh;IScheck et al.ld2004l) ; lBurrows et al.ld2006[|2007l) . and 
IScheck et al. I (120061 henceforth Paper I), a slower development 
of the explosion and the use of a full 180° grid allowed for 
the formation of pronounced global (dipolar and quadrupolar) 
modes of asymmetry. 

The anisotropy in these models is of particular interest, 
as it might provide the explanation f or two results from ob- 
servations: Firstly, s pectropolarimetry (Wang et al. 2001, 2003; 
iLeonard et alj r2006, and references therein) revealed that a non- 
spherical ejecta distribution is a common feature of many core- 
collapse supernovae and is probably caused by the explosion 
mechanism itself, since the anisotropy increases if deeper lay- 
ers of the ejecta are probed. In the case of Supernova 1987 A 
this non-spherical distribution of the ejecta c an even be directl y 
imaged with the Hubble Space Telescope dWang et"aT] 12002). 
Secondly, neutron stars move through interstellar space with 
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velocities much hi gher than those of thei r progenitors (e.g.. 



Cord es et al.l 1 1 993b iLvne & Lorimer 1994; iHansen & Phinney 
1997; Arzoumanian et al. 2002t lZou et alj2005tlChatteriee et all 
2005; iHobbs et al.l 120051) . in some cases with more than 



1000 km/s. It was sugg ested by |Herantl(ll995l) and demonstrated 
with hydrodynamic simulations by Scheck et al.l d2004 [2006) 
that neutron star velocities of this magnitude can result from 
strongly anisotropic (in the most extreme cases "one-sided" i.e., 
dipole-dominated) explosions, in which the total linear momen- 
tum of the ejecta must be balanced by a correspondingly high 
recoil momentum of the neutron star. 

In multi-dimensional simulations convective motions break 
the initial global sphericity and support the explosion (or bring 
the model closer to the explosion threshold) by transporting cool 
matter from the shock to the gain radius where neutrino heating 
is strongest and by allowing hot matter to rise and to increase 
the pressure behind the stalled shock. However, it is not clear 
whether convection can also be responsible for the development 
of low modes in th e postshock accre tion flow, as suggested by 
iHerantj d 19951) and iThompsonl (120001) . The / = 1 pattern stud- 
ied by iHerantl (1 19951) was motivated by a perturbati on analysis 
of volu me-filling convection in a fluid sphere by IChandrasekharl 
(1961), who fo und the dipole (/ = 1) m ode to be the most unsta- 
ble one. In fact. I Woodward et alJd2003l) and lKuhlen etafl d2003h 
demonstrated with three-dimensional simulations that the I = 1 
mode dominates the conv ection in red-giant and main-sequence 
stars. iBlondin et al.l d2003l) . however, investigating an idealized 
setup in 2D hydrodynamic simulations, discovered that an adi- 
abatic accretion flow below a standing shock develops a non- 
radial, oscillatory instability, which they termed "standing accre- 
tion shock instability" or SASI, and which is dominated by the 
I — 1 or / = 2 modes. This suggests that the low-mode asymme- 
tries found to develop in supernova cores in multi-dimensional 
models m ay be caused by glob al instabilities different from con- 
vection. Foglizzo et al. ( 2006) performed a linear stability anal- 
ysis for a problem that resembles the stalled shock situation in 
supernovae, taking into account the limited radial size of the con- 
vectively unstable layer below the shock and the finite advection 
of matter through this region. The latter process turns out to have 
a stabilising effect and can hamper the growth of convection sig- 
nificantly. In particular, the lowest modes are convectively un- 
stable only if the ratio of the convective growth timescale to 
the advection ti me t hrough the unstable layer is small enough. 
Foglizzo et al. (2006) estimate that this may not be the case 
in general and support the suggestion that instabilities different 
from convection may be responsible for the occurrence of low- 
order modes of asymmetry in the postshock accretion flow. 

The "advectiye-acoustic cycle", in short AAC 
(Fog lizzo & Tagger 2000; Foglizzo 2001, 2002), is a promising 
candidate for explaining such a (SASI) instability. It is based on 
the acoustic feedback produced by the advection of entropy and 
vorticity perturbations from the shock to t he forming neutron 
star. B y means of linear stability analysis, iGalletti & Foglizzol 
(2005) showed that due to the AAC the flow in the stalled 
accretion shock phase of core-collapse supernovae is unstable 
with respect to non-radial perturbations, and that the highest 
growth rates are found for the lowest degree modes (in particular 
for the I = 1 mode). 

The situation studied by IBlondin et alj d2003l) and 
IGalletti & Foglizzol d2005l) was, how ever, strongly s impli- 
fied compared to real supernovae. IBlondin et al.l |2003) 
observed the growth of non-radial perturbations in a flow 
between an accretion shock and an inner boundary, which was 
located at a fixed radius. The boundary conditions were taken 



from a stationary flow solution. Furthermore, neither a realistic 
description of the equation of state of t he gas nor the effect s 
of neutrinos were ta ken into account b y Blond irTet alj d2003l) . 
Improving on this, IBlondin & Mezzacappal d2006l) adopted 
an a nalyti c neutrino cooling function ( Houck & Chevalier! 
1992), and lOhnishi et al.l d2006l) in addition took into account 
neutr ino heating and u sed the more realistic equation of state 
from IShen et alj d 1998b . Both groups concur in that low-mode 
instabilities develop also in these more refined simulations. The 
nature of the insta bility mechanism i s, however, still a matter 
of debate. While Ohni shi et alj d2006) consider the AAC as the 
cause of the low-mode oscillations, IBlondin & Mezzacappa 
(2006) argue that a different kind of instability, which is purely 
acoustic and does not involve advection, is at work in their 
simulations. Yet, the eigenmodes found i n the latter simulation s 
were also reproduced in a linear study of Fogliz zo"et alj d2007l) . 
who demonstrated that at least for higher harmonics t he insta- 
bility is the consequence of an advective-acoustic cycle. Laming 
d2007l) . finally, suggested the possibility that feedback processes 
of both kinds can occur and differ in dependence of the ratio of 
the accretion shock radius to the inner boundary of the shocked 
flow near the neu tron star surface. 

Th e wor k by IBlondin & M ezzacappa (2006). lOhnishi et ail 
d2006l) . and iFoglizzo et al.l d2007l) shows that non-radial SASI 
instability of the flow below a standing accretion shock occurs 
also when neutrinos (which could have a damping influence) are 
taken into account. This is in agreement with a linear stability 
analys is of the stationary accretion flow bv lYamasaki & Yamadal 
d2007l) . who included neutrino heating and cooling, and stud- 
ied the influence of varied neutrino luminosities from the proto- 
neutron star. They found that for relatively low neutrino lumi- 
nosities the growth of an oscillatory non-radial instability is fa- 
vored, with the most unstable spherical harmonic mode being a 
function of the luminosity, whereas for sufficiently high neutrino 
luminosity a non-oscillatory instability grows. They attributed 
the former to the AAC and the latter to convection. 

All these studies concentrated on steady-state accretion 
flows, made radical approximations to the employed neutrino 
physics, and considered idealized numerical setups with special 
boundary conditions chosen at the inner and outer radii of the 
considered volume. Because of these simplifications such stud- 
ies are not really able to assess the importance of the different 
kinds of hydrodynamic instabilities for supernova explosions. 
The growth rates of these instabilities depend on the properties 
of the flow, and are thus constant for stationary flows. In real 
supernovae, however, the flow changes continuously, because 
the shock adapts to the varying mass accretion rate, the neutrino 
heating below the shock changes, and the proto-neutron star con- 
tracts. Therefore the growth rates also vary, and a priori it is not 
clear whether they will be high enough for a long enough time 
to allow a growth of some instability to the nonlinear phase on 
a timescale comparable to the explosion timescale (which itself 
can be influenced by the instability and is a priori also unknown). 

The aim of this work is therefore to go some steps further in 
the direction of realism and to abandon the assumption of a sta- 
tionary background flow. To this end we study here the growth 
of hydrodynamic instabilities in a "real" supernova core, i.e., 
we follow in 2D simulations the post-bounce evolution of the 
infalling core of a progenitor star as provided by stellar evo- 
lution calculations, including a physical equation of state for 
the stellar plasma and a more detailed treatment of the neu- 
trino physics than employed in the previous works. The consid- 
ered models were computed through the early phase of collapse 
until shortly after bounce by using state-of-the-art multi-group 
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neutrino transport dBuras et al.ll2003l l2006allbb . In the long-time 
post-bounce simulations performed by us, we then used an ap- 
proximative description of the neutrino transport based on a gray 
(but non-equilibrium) integration of the neutrino number and 
energy equations a long characteristics (for details of the neu- 
trino treatment, see Sch eck et alj |2006 [Paper I]). Compared to 
supernova simulations with a state-of-the-art energy-dependent 
description o f the neutrino t ransp o rt (in spherical 



try see, e.g., [R ampp & Jan kal 120021: [L iebendo rfer et 



1 symme- 
"aiI l200U 



Liebendorfer et al. 2005; Thomp son et al.ll2003l and for multi 
group transport also in 2D. see. e.g jBuras et al.ll2003ll2006albl) 
the models presented here thus still employ significant sim- 
plifications. Such an approximative neutrino treatment must 
therefore be expected to yield results that can differ quanti- 
tatively from those of more sophisticated transport schemes. 
Nevertheless our approach is able to capture the qualita- 
tive features of the better treatments. It is certainly signif- 
icantly more elaborate (and "realis tic") than the sche matic 
neutrino source terms employ ed by IFoglizzo et alj d2007l) and 
iBlondin & Mez zacappa (2006), and the l ocal neutrino source 
desc ription (without transport) a dopted bv lOhnishi et al.l (2006) 
and lYamasaki & Yamadal d2007l) . We consider our approxima- 
tion as good enough for a project that does not intend to establish 
the viability of the neutrino-heating mechanism but which is in- 
terested mostly in studying fundamental aspects of the growth of 
non-radial hydrodynamic instabilities in the environment of su- 
pernova cores including the influence that neutrino cooling and 
heating have in this context. 

We made use of one more approximation that reduces the 
complexity of our simulations compared to full-scale supernova 
models, namely, we did not include the neutron star core but re- 
placed it by a Lagrangian (i.e., comoving with the matter) inner 
grid boundary that contracts with time to smaller radii, mim- 
icking the shrinking of the cooling nascent neutron star. At this 
moving boundary the neutrino luminosities produced by the neu- 
tron star core were imposed as boundary conditions. This had 
the advantage that we could regulate the readiness of a model 
to explode or not explode, depending on the size of the cho- 
sen core luminosities and the speed of the boundary contrac- 
tion. The inner boundary of our computational grid is impen- 
etrable for the infalling accretion flow, but the accreted matter 
settles into the surface layer of the forming neutron star, sim- 
ilar to what happens outside of the rigid core of the compact 
remnant at the center of a supernova explosion. This is different 
from the various kin ds of "outflow bound aries " employed in the 
literatureQ although IBlondin et al.l d2003l) and lBlondin & Shawl 
(2007) reported about tests with several different prescriptions 
for the boundary treatment without finding any significant influ- 
ence on the growth of the SASI. Our modeling approach there- 
fore follows Schec ket al.l (120041) and Paper I, where indeed the 
development of low-mode flow (with dominant / = 1 and I = 2 
modes) between shock and neutron star was found. In these pre- 
vious papers we, however, did not attempt to identify the mech- 
anism^) that were causal for the observed phenomenon and just 
mentioned that convection and the acoustically-driven or AAC- 
driven SASI may yield an explanation for the large global asym- 



1 IBlondin et alj j2003l) used a "leaky boundary" and lOhnishi et alj 
d2006h a ' 'free outflow boundary", bot h assuming non-zero radial ve - 
locity at the grid boundary. In contrast, Blondin & Mezzacappa ( 2006) 
adopted a "hard reflecting boundary". Although in this case the radial 
velocity at the boundary is taken to be zero, there is still a non-zero 
mass flux as the density near the boundary goes to infinity. The differ- 
ence can be important for the spurious generation of acoustic feedback 
by vorticity perturbations. 



metries seen to develop during the neutrino-heating phase of the 
stalled shock. There was no analysis which mechanism was ac- 
tive and why it had favorable conditions for growth. 

In the present work we return to these questions. In particular 
we aim here at exploring the following points: 

- What is the timescale for a non-radial instability of the 
stalled accretion shock (SASI) to develop in a supernova 
core? How is it influenced by neutrino effects? 

- Can the instability be identified as consequence of an ampli- 
fying adve ctive-acoustic cycle, of a gro wing standing pres- 
sure wave dBlondin & M ezzacappa 20061) or of something 
else? 

- What determines its growth rate? For which conditions does 
the instability grow faster than convection and which influ- 
ence may this have on the subsequent (nonlinear) evolution? 

- What is the relationship between convection and the insta- 
bility in the nonlinear phase? 

- What is the possible supportive role of the SASI in the con- 
text of neutrino-driven explosions and in creating the low- 
mode ejecta asymmetry identifie d as cause of large neutr on 
star kicks bv lScheck et ail d2004 and lScheck etall d2006l) ? 

In order to address these questions we will first summarise 
the most important properties of convection and of the AAC in 
the gain layer in Sect. [2] This will serve us as basis for the later 
analysis of our set of two-dimensional simulations. In Sect.[3]we 
will describe the computational methods and numerical setup we 
use for these simulations and will motivate our choice of param- 
eter values for the considered models. We will present the simu- 
lation results in Sect.|4]and will discuss them in detail in Sects. [5] 
and[6] Section[7] finally, contains our conclusions. 

2. Hydrodynamic instabilities 

2. 1 . Linear and nonlinear convective growth of perturbations 

In a hydrostatic, inviscid atmosphere, regions with negative en- 
tropy gradients (disregarding possible effects of composition 
gradients) are convectively unstable for all wavelengths. Short 
wavelength perturbations grow fastest, with a local growth rate 
i^buoy > equal to the imaginary part of the complex Brunt- 
Vaisala frequency: 



k^buoy 



'grav 



C/ P . 



(1) 



Here a glav < is the local gravitational acceleration, p is the 
density and 



as 



P,Y, 



dr 



(2) 



where S is the entropy, P is the pressure, and F/ is the total lepton 
number per nucleon. Note that C > is the instability condition 
for Schwarzschild convec tion. 

IFoglizzo et al.l (|2006) pointed out that in the stalled shock 
phase, the convective growth timescale 0J b ^ 0y in the unstable 
layer below the shock is of the same order as the timescale for 
advection from the shock to the gain radius, 



•r 



dr 



(3) 



where R g is the gain radius, R s the shock radius and v r the ra- 
dial velocity. Advection is stabilising because it gives pertur- 
bations only a finite time to grow in the gain region, before 
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they are advected into the stable layer below the gain radius. 
Considering the local growth rate a>buoy(/) given by Eq. ([TJ in 
a reference frame advected with the flow, the amplitude 8 of 
a small-wavelength perturbation may grow during its advection 
from the shock to the gain radius, at best by a factor exp(y), 



<5gain = £ s hock ex P(^)> 



where the quantity 

X = I w buoy (r) 



dr _ g 
v r (r)| ~ T ^' T ° 



(4) 



(5) 



can be interpreted as the ratio of the advection timescale to 
the average local growth timescale the perturbation experiences, 
i"conv = (^buoy) ( tne latter quantity is implicitly defined by Eq.|5]). 
Thus it would appear that in order to reach a given perturbation 
amplitude at the gain radius, a certain seed perturbation ampli- 
tude of the matter crossing the shock would be necessary. 

However, a linear stability analysis reveals that the station- 
ary accretion flow below the shock is globally unstable and 
perturbations can grow from arbitrarily small in itial seeds, if 
suffic ient time is availa ble (Foglizzo et al.l [2006b . According 
to iFoglizzo et a D 42006) this is the case for a limited range 
['min, 'max] of modes for which ^ exceeds a critical value xo, 



X>Xo 



where 



Xo 



(6) 



For^- < xo the flow remains linearly stable, even though a nega- 
tive entropy gradien t is present. 

The analysis of IFoglizzo et al.l (120061) applies only for the 
linear phase of the instability, i.e. for small perturbation ampli- 
tudes. However, it is possible that the situation has to be con- 
sidered as nonlinear right from the beginning, i.e. that the seed 
perturbations grow to large amplitudes already during their ad- 
vection to the gain radius. In this context "large" can be defined 
by considering the buoyant acceleration of the perturbations. 

For a small bubble, in which the density p is lower than the 
one of the surrounding medium, p surr , the convective growth dur- 
ing the advection to the gain radius may lead to an increase of the 
relative density deviation 5 =\p - p SU n-|/PsuiT (which can be con- 
sidered as the perturbation amplitude) as given by Eq. ©. The 
bubble experiences a buoyant acceleration |a gla vl<5 towards the 
shock, which is proportional to the local gravitational accelera- 
tion flgrav The time integral of the buoyant acceleration becomes 
comparable to the advection velocity, when the perturbation am- 
plitude reaches a critical value 

(Ifrl)gain 

Ocrit = ; 5— (7) 



<a g rav> 



grav/garn ' adv 



<|Vrl>: 



gain 



(^grav)gain ^g 



0(1%), 



(8) 



where (|v r |) gam and (a grav )gain are the average values of the radial 
velocity and the gravitational acceleration in the gain layer, re- 
spectively. For (5gaj n > (J cr it a small-scale perturbation is able to 
rise against the accretion flow. If the whole flow is perturbed, the 
buoyant motions on small scales affect the situation globally and 
could allow for the onset of convective overturn also on larger 
scales. Note that in contrast to the linear growth of the instabil- 
ity this process does not require ;f > xo but it does require large 
enough seed perturbations, 

c ^crit /r .* 

expCr) 

A sufficient condition for the suppression of convection is there- 
fore that neither Eq. (|6]l nor Eq. (O are fulfilled. 
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Fig. 1. Schematic view of the advective-acoustic cycle between 
the shock at R s (thick solid line) and the coupling radius, R c 
(thick dashed line), in the linear regime, shown for the case 
where the oscillation period of the shock (t osc ) equals the cycle 
duration, r aac . Flow lines carrying vorticity perturbations down- 
wards are drawn as solid lines, and the pressure feedback cor- 
responds to dotted lines with arrows. In the gray shaded area 
around R c the flow is decelerated strongly. 



2.2. The advective-acoustic cycle 

A second hydrodynamic instability has recently been recognised 
to be of potential importance in the stalled shock phase. Blondin 
et al. (2003) noticed that the stalled accretion shock becomes un- 
stable to non-radial deformations even in the absence of entropy 
gradients, a phenomenon termed SASI. It can be interpreted as 
the result of an " advective-acoustic cycle" ( in short AAC), as 
first discussed by Foglizzo & Tagger (2000) in the context of 
accretion onto black holes , and later studied for sup ernovae by 
iGalletti &Foglizzol d2005l) and IFoglizzo et al.l J2007) by means 
of linear stability analysis. The explanation of these oscillations 
is based on the linear coupling between advected and acoustic 
perturbations due to flow gradients. 

Although this linear coupling occurs continuously through- 
out the accretion flow from the shock to the neutron star surface, 
some regions may contribute more than others to produce a pres- 
sure feedback towards the shock and establish a global feedback 
loop. The analysis of the linear phase of the instability in Sect. [5] 
reveals the importance of a small region at a radius R c above the 
neutron star surface, where the flow is strongly decelerated. The 
feedback loop can be described schematically as follows: small 
perturbations of the supernova shock cause entropy and vorticity 
fluctuations, which are advected downwards. When the flow is 
decelerated and compressed above the neutron star surface, the 
advected perturbations trigger a pressure feedback. This pressure 
feedback perturbs the shock, causing new vorticity and entropy 
perturbations. Instability corresponds to the amplification of per- 
turbations by a factor I2aacl > 1 through each cycle. 

The duration r aac of each cycle is a fundamental timescale. 
It corresponds to the time needed for the advection of vortical 
perturbations from the shock to the coupling radius R c , where 
the pressure feedback is generated, plus the time required by the 
pressure feedback to travel from this region back to the shock. 

The oscillatory exponential growth resulting from the AAC 
can be described by a complex eigenfrequency co = co r + iu>j 
satisfying the following equation: 



expt-zW^) = a 



(10) 
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A fl/f 
AM gain 




A/t 
«11S 






[B/s] 


[km] 


[s] 




[B/s] 


[s] 


[s] 




[B] 


[M e ] 


[km/s] 


wnnin a 
wuurA 




O.U 


U.J 


1 n~ 3 1, 


u.u 














WOOF 


0.2 


8.0 


0.5 


W~ 3 v r 


99.0 


0.154 


0.194 


0.004 


0.37 


1.50 


200 


WOO 


0.2 


15.0 


1.0 


\0- 3 V r 


51.7 


0.346 








1.63 


-3 


woos 


0.2 


15.0 


4.0 


10" 3 V,- 


29.6 














W05S 


7.4 


15.0 


4.0 


10 6 cm/s 


38.3 














W05V 


7.4 


15.0 


10.0 


3 x 10 7 cm/s 


33.5 














W12F 


29.7 


10.5 


0.25 


10- j v r 


109.3 


0.144 


0.164 


0.010 


0.87 


1.44 


-558 


W12F-C 


29.7 


10.5 


0.25 


0(lQ- 2 )v r 


112.5 


0.090 


0.117 


0.015 


0.94 


1.41 


612 



' The constant boundary luminosity is denoted by Lib, R\ h is the asymptotic inner boundary radius (see Eq.ll3t. fib the contraction timescale of the 
inner boundary, 6i the amplitude of the initial velocity perturbations (which was chosen to be proportional to the local radial velocity v, in most 
cases — see text for details), L Vl . +r . L . is the luminosity of electron neutrinos and antineutrinos at a radius of 500 km at 150 ms after the start of the 
simulations, t nl the time at which the average lateral velocity in the gain layer exceeds 10 9 cm/s, AM gain the mass in the gain layer at the time ? exp 
when the explosion starts, E ew the explosion energy, M ns the neutron star mass, and v ns the neutron star velocity as computed from momentum 
conservation. The last three quantities are given at t = Is, for Model WOOF at t = 750 ms, and / cxp is denned as the moment when the energy of the 
matter in the gain layer with positive specific energy exceeds 10 49 erg. Model W00FA is a hydrodynamic simulation without including neutrino 
effects. Only Models WOOF, W12F, and W12F-C developed explosions and values for the corresponding explosion and neutron star parameters are 
given. The neutrino luminosities imposed at the inner boundary are kept constant during a time t L = 1 s. The energy unit 1 B = 10 51 erg is used 
and all times are measured from the start of the simulation, i.e. t = s means 16 ms after core bounce. 



where the real part oj r is the oscillation frequency and the imag- 
inary part cji is the growth rate of the A AC. Not e that Eq. ( fTQb is 
a simplified form of Eq. (49) of Foglizzo] d2002l) . For the sake of 
simplicity, it neglects the marginal influence of the purely acous- 
tic cycle of pressure waves trapped between the shock and the 
accretor. This hypothesis is motivated by the recent estimate of 
the efficiencies Q mc , 7? ac of the advective-acou stic and of the 
purely acoustic cycles, respectively, obtained by Foglizzo et al. 
d2007l) in a simpler set-up when the frequency is high enough to 
allow for a WKB approximation. According to their Figs. 8 and 
9, the purely acoustic cycle is always stable (|^ ac | < 1), whereas 
the advective-acoustic cycle is unstable (IQaad > lfl According 
to Eq. ( [Tol l, the oscillation period t osc = 2n/a> r of the AAC de- 
pends both on the duration T aac of the cycle and on the phase <p 

of gaa C : 



The amplitude of perturbations in the AAC increases like 
exp(w, f), with a growth rate <w, deduced from Eq. ( fTOl : 



w r"r aac + <p = 2nn, 



(11) 



where n is an integer lab elling the diffe r ent ha rmonics. In the 
particular flow studied by Fogliz zo"et al.l d2007l) . the oscillation 
period of the fundamental mode is a good estimate of the du- 
ration of the cycle (t osc ~ T aac ). This simple relationship is not 
obvious a priori. For example, if (2 aac were real and negative (i.e., 
ip = 7r), the oscillation period would scale like r osc ~ 2r aac (be- 
cause enhancing feedback requires a phase coherence between 
the amplifying mechanism and the shock oscillation). 



2 This discussion applies for the growth behavior of the advective- 
acoustic cycle and of the purely acoustic cycle in the presence of a pres- 
sure feedback produced at the coupli ng radius and reaching the shock. 
The calculation of 7? ac and Qaac by iFoglizzo et alj d2007l) does not as- 
sume a purely radial acoustic feedback but fully takes into account the 
azimuthally traveling sound waves as well as evanescent pressure waves 
(pseudosound), which do not propagate. Ne vertheless, it is c urre ntly a 
controversial issue whether the analysis by Foglizzo et al.1 d2007l) al- 
lows one to draw concl usions on th e kind of instability p roposed by 
Blondin & Mezzacappa ( 2006) and Blondin & Shaw ( 2007), who advo- 
cate a growth mechanism driven by sound waves traveling solely in the 
angular direction. The possibility of understanding the development of 
SASI modes in our hydrodynamic supernova simulations by an acous- 
tic cycle with non-radial sound wave propagation will be discussed in 
Sect.[U 



lndQaacI) 



(12) 



Comparing Eqs. (O and (fT2l it is interesting to note that a small 
advection timescale suppresses the growth of entropy-driven 
convection whereas it leads to higher growth rates for the AAC 
(neglecting the logarithmic dependence on £) aac ). Thus the AAC 
may operate under conditions which are not favourable for con- 
vection and vice versa. Investigating this further is one of the 
goals of this work. 



3. Numerical setup and models 

In order to investigate the importance of instabilities like the 
ones discussed in the previous section during the post-bounce 
evolution of core-collapse supernovae, we performed a series 
of two-dimensional (2D) hydrodynamic simulations. For this 
purpose we used the same numerical setup as in Paper I. We 
employed the version of the hyd rodynamics code that was de- 
scribed by Kifonid is et al.l d2003l). It is based on the piec ewise 
parabolic method (PPM) of IColella & Woodward! (1 19841) , as- 
suming axisymmetry and adopting spherical coordinates (r, 9). 
The calculations were performed on a polar grid that had typi- 
cally 800 zones in radial direction and 360 zones in lateral di- 
rection (extending from polar angle 8 = to 6 — n). The lat- 
eral grid was equidistant while the radial grid had logarithmic 
spacing with a ratio of radial zone size to radius that did not 
exceed 1%. For the neutrino number and energy transport we 
applied a gray, characteristics-based transport scheme that was 
able to efficiently approximate the transport in the transparent 
and semi-transparent regimes up to optical depths of several 100. 
Only transport in the radial direction was taken into account, but 
we allowed for lateral variations of the neutrino flux by solving 
one-dimensional transport equations independently for all dis- 
crete polar angles of the r-9 grid. A detailed description of the 
transport method is given in Paper I. 
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Fig. 2. Entropy distribution of model WOOFA 30 ms and 190 ms 
after the start of the simulation. The initial entropy profile and 
postshock entropy gradients caused by shock motions give rise 
to weak convection. A low-amplitude / = 1 oscillation develops. 
( Color figures are available in the online version of our paper.) 
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Fig. 3. Mass-shell trajectories for model WOOFA. The spacing 
of the thin lines is 0.01 M Q . Green lines mark the mass shells at 
which the composition of the progenitor changes. The red lines 
are the minimum, average, and maximum shock radii, the black 
line marks the radius, at which the average density is 10 11 g/cm 3 . 
The difference between minimum and maximum shock radius is 
caused by bipolar shock oscillations (see Fig. 13. (A color figure 
is available in the online version of our paper. ) 



3. 1 . Initial and boundary conditions 

We used an initial model (the 'W' model from Paper I) that was 
o btained by evolving t he 15 M supernova progenitor sl5s7b2 
of IWooslev & Weaver! £f995) through collapse and core bounce 
until shock stagnation in a simulation wit h a detailed, en ergy- 
dependent treatment of neutrino transport dBuras et al. 20031 see 



their Model si 5). We started our runs at a time of 16 ms after 
core bounce. In order to enable the growth of hydrodynamic in- 
stabilities we perturbed the initial model, unless noted otherwise, 
by adding random, zone-to-zone velocity perturbations of 0. 1 % 
amplitude. 

The neutron star core (i.e. typically the innermost 1 . 1 M ) 
was not included in our simulations but was replaced by a con- 
tracting inner boundary of the computational grid. Boundary 
conditions were imposed there for the hydrodynamics and the 
neutrino transport, and a point-mass potential of the excised core 
was adopted to account for the gravitational influence of this re- 
gion. Although the treatment of gravity is not of primary rele- 
vance for the fundamental questions studied in this paper, we 
mention here that the description of the gravitational potential 
took into account the self-gravity of the gas on the grid with its 
two-dimensional distribution, as well as an approximative treat- 
ment of general relativistic effects (for details, see Paper I). 

The inner grid boundary was placed at a Lagrangian shell 
with enclosed mass of M = 1.1 M at which we imposed condi- 
tions describing hydrostatic equilibrium. Its radius was assumed 
to evolve according to 

R- h 

Rib® = j f , (13) 

1 + (l-exp(-f/f ib ))(^ b /^ b -D 

where Rt is the final (asymptotic) boundary radius, t^ is the 
contraction timescale and R\ h « 65 km is the initial radius, which 
is given by the initial model. The neutrino luminosities from the 
neutron star core, which we imposed at the inner boundary, were 
assumed to be constant during the first second after core bounce. 
This simple choice can be justified by the results of core-collapse 
simulations with sophisticated neutrino transport (see Paper I). 

With this approach we parametrized the cooling and shrink- 
ing of the core of the nascent neutron star and its neutrino emis- 
sion, which all depend on the incompletely known properties of 
the nuclear equation of state. Different choices of the bound- 
ary motion and strength of the neutrino emission allowed us to 
vary the properties of the supernova explosion and of the devel- 
oping hydrodynamic instabilities in the region between neutron 
star and stalled shock. It is very important to note that the stag- 
nation radius of the stalled shock reacts sensitively not only to 
the mass infall rate from the collapsing progenitor star and to the 
rate of neutrino heating in the gain layer, but also to the contrac- 
tion behavior of the neutron star. A faster contraction usually 
leads to a retraction of the shock, whereas a less rapid shrink- 
ing of the neutron star allows the shock to expand and stagnate 
at a larger radius. This, of course, causes important differences 
of the postshock flow and thus affects the growth of non-radial 
hydrodynamic instabilities. 

In some of the simulations discussed here, the rapid contrac- 
tion of the forming neutron star caused the density and sound 
speed at the inner boundary to become so high that the hydrody- 
namic timestep was severely limited by the Courant-Friedrich- 
Lewy (CFL) condition. Moreover, when the optical depth in this 
region increased to more than several hundred, numerical prob- 
lems with our neutrino transport method occurred unless very 
fine radial zoning was chosen, making the timestep even smaller. 
In such cases we moved the inner grid boundary to a larger ra- 
dius and bigger enclosed mass (i.e., we increased the excised 
neutron star core). Hereby we attempted to change the contrac- 
tion behavior of the nascent neutron star as little as possible. The 
new inner boundary was placed at a radius Rfa where the opti- 
cal depth for electron neutrinos was typically around 100. When 
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doing this, the gravity-pr oducing mass of the inner core was ad- 
justed appropriately (see lArcones et al. 2006) and the boundary 
neutrino luminosities were set to the values present at R{\, at the 
time of the boundary shifting, thus making sure that the gravi- 
tational acceleration, the neutrino flux, and neutrino heating and 
cooling above R(\, followed a continuous evolution. The parame- 
ters in Eq. ( TT~3T > were adjusted from the old values R 

new values RL and /?£, respectively, in the following way: 



ib and R jh to 



4, 



R 



R^xiRURib), 
I = Rl x (R( b /R ib ) ■ 



(14) 



This simple rescaling had the consequence that during the sub- 
sequent evolution small differences of the contraction velocity 
of the new inner boundary and therefore of the settling neutron 
star appeared, which led also to minor changes of the decay of 
the neutrino luminosities with time. Nevertheless, no significant 
impact on the simulations was observed, e.g., Model WOOF de- 
veloped an explosion at the same time, independent of whether 
or not the boundary was shifted according to our recipe. 
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Fig. 4. Evolution of the quantity^ (upper panel, see Eq.|5j and of 
the ratio <5 ga in /<5crit (lower panel, see Eqs.|4j[8]l for Models W12F- 
c, W12F, WOOF and WOO. All lines end 10 ms before the non- 
linear phase begins at time t n \, when the average lateral velocity 
in the gain layer exceeds 10 9 cm/ s. At later times x an d ^crit can- 
not be measured reliably any longer. In all models, x 5 Xo ~ 3 
at t < t n \. Only in Model W12F-C the ratio <5gain/£crit gets very 
close to unity before t = t n \, which means that only in this model 
convection is able to set in faster than the SASI. 
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Fig. 5. Same as Fig. [3] but for Models WOO (upper panel) and 
WOOF (lower panel). The blue line marks the position of the gain 
radius. Up to f » 300 ms Model WOO remains nearly spherical 
and evolves like the corresponding one-dimensional simulation. 
However, an initially very weak / = 1 oscillation mode in the 
postshock flow grows in this phase and finally becomes non- 
linear, causing strong shock oscillations. Yet, this model does 
not explode. Although the shock expands transiently in a quasi- 
periodic manner, the average shock radius decreases and all mat- 
ter remains bound. In Model WOOF an / = 2 mode develops and 
starts to affect the shape of the shock at t * 150 ms, much ear- 
lier than in Model WOO. The oscillations become nonlinear, and 
at t = 194 ms (marked by a dotted vertical line) the model ex- 
plodes. (Color figures are available in the online version of our 
paper. ) 



°- 4 3.2. Model parameters 



The characteristic parameters and some important quantities of 
the eight models investigated here are listed in TableQ] The mod- 
els differ concerning the included physics, assumed boundary 
conditions, and the initial perturbations used to seed the growth 
of hydrodynamic instabilities. 

The most simplified case we considered, Model W00FA, is 
a purely hydrodynamic simu lation without incl uding neutrino 
effects. This choice follows iBlondin et al] (2003), who also ig- 
nored neutrinos. In comparison with the other models we com- 
puted, it allows us to study th e influence of neutrino cooling and 
heating. IBlondin et alj (12003) also placed an inner boundary at a 
fixed radius and applied outflow conditions there to allow for a 
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steady-state accretion flow (alternatively, they also tested reflect- 
ing conditions with a cooling term to keep the shock at a steady 
radius). In contrast, in Model WOOFA accretion is enabled by the 
retraction of the inner boundary of the computational grid, which 
mimics the Lagrangian motion of a mass shell in a contr acting 
neutron star. Another difference from[Blondin et al. (2003) is the 
fact that in our models the accretion rate shrinks when infalling 
matter from the less dense layers at increasingly larger radii 
reaches the shock. Thus the development of hydrodynamic in- 
stabilities occurs in a situation that is genetically non-stationary. 

In five other simulations we included neutrinos and chose 
boundary conditions such that the growth of convection was sup- 
pressed. This allowed us to identify and study other instabilities 
like the SASI more easily. The suppression of convection could 
be achieved by prescribing vanishing or negligibly low core lu- 
minosities. In such cases only the luminosity produced between 
the inner boundary and the gain radius causes neutrino energy 
deposition in the gain layer. Therefore the neutrino heating re- 
mains weak, resulting in a shallow entropy gradient and con- 
sequently in a large growth timescale for convection. This im- 
plies that for low core luminosities the ratio of the advection 
to the buoyancy timescale, y (Eq. H), remains below the criti- 
cal value and therefore in spite of a negative entropy gradient 
in the neutrino heating region, the postshock layer remains con- 
vectively stable due to the rapid advection of the gas down to 
the gain radius (see Sect. 12. It . The five simulations where this 
is the case are Models WOOF, WOO, W00S, W05S, and W05V. 
These models differ in the prescribed contraction of the inner 
boundary. Models WOO and WOOF employ the "standard" and 
"rapid" boundary contraction, respectively, of Paper I. In order 
to cover a wider range of advection timescales — which will 
help us to gain deeper insight into the mechanism that causes the 
low-mode instability found in our simulations (see Sect.[5} — we 
performed three simulations with slower boundary contraction, 
namely Models W00S, W05S, and W05V (Table Q). In the last 
two models the core neutrino luminosity has a non-negligible 
(but still fairly low) value. The correspondingly enhanced neu- 
trino heating leads to larger shock radii and thus longer advec- 
tion timescales. Models WOOF, WOO, and WOOS were computed 
with our standard initial perturbations (0. 1 % random noise on 
the velocity). For Models W05S and W05V an I = 1 velocity 
perturbation was applied. This allowed us to suppress high-mode 
noise and to measure the oscillation period of the low-mode in- 
stability despite the low growth rates in these models. 

For Models W12F and W12F-C, finally, we adopted bound- 
ary conditions that were guided by core-collapse simulations 
with sophisticated multi-group neutrino transport. The contrac- 
tion of the inner boundary was chosen to match the motion of 
the corresp onding mass shell in such simulations for the same 
progenitor dBuras et al.l l2003). The boundary luminosity we im- 
posed led to typical explosion energies of about 10 51 erg. Despite 
the non-negligible core luminosity convection in these models 
was suppressed because of the rapid boundary contraction. The 
latter caused the radius of the stalled shock to become rather 
small, and consequently the accretion velocities in the postshock 
layer were very large. Therefore the advection timescale was 
short and the parameter x did not exceed the critical value of 
about 3. In such a situation the amplitude of the progenitor per- 
turbations can decide about whether convection sets in (start- 
ing in the nonlinear regime as discussed in Sect. 12. U or not. 
Since the properties of the perturbations in the progenitor star are 
not well known, we decided to explore two cases, one (Model 
W12F) with small initial perturbations (our standard 0.1% ve- 
locity perturbation) such that the growth of convection was sup- 



pressed, and another case where the initial perturbations were 
large enough so that convection could develop. For the latter 
model, W12F-C, we used the same perturbations as for Models 
W12-C and W18-C of Paper I with amplitudes of up to several 
percent and a spatial variation as given by the velocity fluctu- 
ations that had grown during a 2D core-collap se si mulation of 
a 15 M H star (Model sl5r of iBuras et alll2003l and iBuras et all 
20063)- 



4. Results 

In this section we will give an overview of the simulation results, 
whose interpretation will be given in more detail in Sects. [5]and 

m 

4.1. A model without neutrinos 

Although Model WOOFA does not include neutrino heating, con- 
vective fluid motions develop in this case because a convectively 
unstable region with a negative entropy gradient is present at 
r « 70 km already in the initial conditions of our simulations. 
This feature is a consequence of the decreasing shock strength 
before shock stagnation. Soon after we start model model run 
WOOFA, buoyant bubbles form in the unstable region and rise to- 
wards the shock (Fig. |2j. Convective action continues during the 
whole simulation because neutrino cooling, which could damp 
convection, is disregarded, and because convectively unstable 
entropy gradients are created by shock motions that are caused 
by variations of the preshock accretion rate and by bipolar shock 
oscillations due to SASI modes (see Sect. I6.ll ). However, with- 
out neutrino heating the convective overturn does not become 
as strong and dynamical as in the simulations of Paper I, where 
neutrino effects were included. Also the bipolar shock oscilla- 
tions are rather weak (the shock deformation amplitude does not 
exceed 15%) and occur quasi-periodically with a period of 20- 
50 ms (Fig. 0). 

These multi-dimensional processes do not affect the over- 
all evolution of the model and the shock position as a function 
of time is almost identical to the one found in a corresponding 
one-dimensional simulation. In spite of the contraction of the in- 
ner grid boundary, the shock expands slowly and continuously 
(Fig.|3j. A transient faster expansion occurs at t « 150 ms, when 
a composition interface of the progenitor star falls through the 
shock and the mass accretion rate drops abruptly. After 660 ms 
we stopped the simulation. At this time the shock had reached a 
radius of 400 km. 

Although the shock expands slowly, this does not lead to an 
explosion because without neutrino heating the specific energy 
of the matter behind the shock remains negative. The shock ex- 
pansion takes place because matter piles up in the postshock re- 
gion and forms an extended atmosphere around the neutron star. 
This slowly pushes the shock further out in response to the ad- 
justment of hydrostatic equilibrium by the accumulation of mass 
in the downstream region. Since in the absence of cooling pro- 
cesses the matter cannot lose its entropy, it is not able to set- 
tle down onto the neutron star quickly. Therefore the behavior 
of Model WOOFA is destinctively different from the situation 
obtained in supernova simulations with neutrino transport, and 
it also differs from the stationary flow that was considered by 
iBlondin et ail d2003l) . The postshock velocity in Model WOOFA 
is much lower and the shock radius becomes larger. 

Neutrino cooling is therefore essential to obtain a quasi- 
steady state accretion flow when simulations are performed in 



Scheck et al.: Multidimensional supernova simulations 



9 



1.2 

,— , 1-0 
o 0.8 
O 0.6 

1—1 0.4 

x 

0.2 



1.2 
i—i 1.0 

E 

o 0.8 
O 0.6 




-1.0 -0.5 0.0 0.5 
z [10 7 cm] 



1.0 



0.4 
0.2 



1 .4 
,— , 1-2 

o 1-° 

•b °' 8 
0.6 

x 0.4 

0.2 



-0.5 0.0 0.5 
z [10 7 cm| 



; ; ;t; <=' t v',z,mWtC 



-0.5 0.0 0.5 
z [10 7 cm] 



I 





1.0 
0.8 
0.6 
0.4 
0.2 




-0.5 0.0 0.5 

z [10 7 cm] 



E 

o 

ID 
O 




-5 

z [10 6 cm] 



1.4 
-, 1.2 
o '-01 

0.8 
Zj 0.6 
x 0.4 

0.2 



^OO-P" 'lateral 1 velocity pvd 9 ' crn/ / s ; F- 

\\\\\\\\ \ \ \ \ i ; i j i i f J V / / / 7 S S S 4* 




z [10 7 cm] 



Fig. 6. Lateral velocities (color coded; superimposed are the vectors of the velocity field, which indicate the direction of the flow) 
for Models WOO and WOOF. The white lines mark the shock, the black dotted lines the gain radius. For both models we show the 
situation at two times near t = t n \, at which the oscillations are in opposite phases (i.e. the times differ by half an oscillation period), 
and at a third time, when the oscillations have run out of phase (see Sect. 16. II ). In Model WOO an I — 1 SASI mode develops, i.e., the 
still nearly spherical shock moves back and forth along the z-axis, whereas in Model WOOF an I — 2 SASI mode becomes dominant, 
i.e., the shock oscillates between a prolate and an oblate deformation. The postshock matter attains high lateral velocities, because 
the radial preshock flow hits the shock at an oblique angle when the shock is nonspherical or when its center is displaced from the 
grid center. ( Color figures are available in the online version of our paper. ) 



which the central neutron star is included (in our models it is 
partly excised and replaced by an impenetrable inner grid bound- 
ary). Only when neutrinos remove energy and reduce the entropy 
of the gas can the matter be integrated into the dense surface 
layers of the compact object. The rapid flow of the gas from 
the shock to the neutron star implies short advection timescales, 
which are crucial for the growth of the SASI (see the discus- 
sion in Sec. 12.21 ). Although the accretion flow that develops 
in our supern o va si mulation s is similar t o the o ne assumed by 
iBlondin et aT] d2003l) and lOhnishi et all d2006l) . there are still 
potentially important differences. Because of the contraction of 
the neutron star and due to the density gradient in the collaps- 
ing star, the mass accretion rate varies (usually decreases) with 
time and the accretion between shock and neutron star surface 
never becomes perfectly stationary . Ou r simulations also differ 
from t hose of lOhnishi et al.l d2006l) and lBlondin & Mezzac appa 
(2006) by our more detailed treatment of the neutrino effects. 
Altogether this allows us to assess the questions how non-radial 
hydrodynamic instabilities develop at more realistic model con- 
ditions for the supernova core than considered in previous stud- 
ies, and how such instabilities may influence the onset of the 
supernova explosion. 



4.2. Models with suppressed convection 

In the models including neutrino transport the accreted matter 
loses energy and entropy by neutrino cooling and thus is able to 
settle down onto the neutron star, following the contraction of 
the inner boundary. Comparing the mass shell trajectories of the 
neutrinoless Model W00FA and of Model WOO (shown in Figs. [3] 
and |5] respectively) this difference becomes evident. Since the 
accreted matter does not pile up, also the shock turns around af- 
ter an initial expansion phase and recedes continuously during 
the later evolution (except for a short, transient expansion phase 
a t a 150 ms, which is initiated when a composition interface 
of the progenitor star crosses the shock). Due to the miniscule 
boundary luminosity the neutrino heating remains weak and the 
parameter x °f Eq. © stays below the critical value (Fig. |4]). 
Consequently, there is no evidence of convection in the gain 
layer and Model WOO evolves nearly spherically symmetrically 
in the first 300 ms. 

However, already several ten milliseconds after the start of 
the simulation a lateral velocity component (which changes di- 
rection with a period of about 30 ms) is observable in the flow be- 
tween shock and neutron star surface. The amplitude of this I = 1 
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Fig. 7. Entropy distribution of Model WOO for several moments near the beginning of the nonlinear phase (the displayed times have 
a separation of half an oscillation period), and at t — Is. Within each SASI oscillation cycle the postshock entropies vary strongly 
and steep, unstable entropy gradients develop in the postshock flow. Finally, the Rayleigh-Taylor growth timescale becomes smaller 
than the oscillation period and the characteristic mushroom structures are able to grow. In the subsequent evolution the low-mode 
oscillations saturate and the model does not develop an explosion. ( Color figures are available in the online version of our paper. ) 



oscillation mode starts to increase continuously after t » 100 ms 
and grows by a factor of about two per period. However, the 
amplitude is not large enough to affect the shape of the shock 
before t « 250 ms because the finite resolution of the numeri- 
cal grid prevents the shock from being pushed out by less than 
one radial zone and thus it remains perfectly spherical for low 
oscillation amplitudes (lateral variation is already visible in the 
postshock flow, though). 

In the subsequent evolution the shock radius is initially still 
slowly decreasing and the shock shape remains approximately 
spherical, but the shock surface moves back and forth along the 
axis of symmetry assumed in our two-dimensional simulations. 
The direction of the postshock flow changes periodically and 
the flow transports matter between the southern and the northern 
hemispheres (Fig. |5J. This situation is quite similar to the bipo- 
lar oscillations encountered in some models discussed in Paper I 
and also in full-sca le supernova simul ations with sophisticated 
neutrino transport (Bur as et"aill2006bj) . However, because con- 
vection is absent, the flow pattern and the shape of the shock are 
much less structured in Model WOO. 

At t « 360 ms the amplitude of the shock oscillations has 
become very large, the shock radii at the poles differ by up to 
50 km, whereas the average shock radius is only about 100 km. 
In this phase the entropy behind the shock starts to vary strongly 
with time and angle (Fig. Q. Steep negative entropy gradients 
{AS /dr — 0(1 £t,/km)) develop and Rayleigh-Taylor instabilities 
start to grow at the boundaries between low- and high-entropy 
matter. The postshock flow reaches lateral velocities of several 
10 9 cm/s and supersonic downflows towards the neutron star 
form (see Sect. 16. II for a discussion of these processes). Within 
a few oscillation cycles the whole postshock flow becomes very 
similar to the nonlinear convective overturn present at the onset 
of the explosion in those models of Paper I where the explosion 
energy was rather low. 

However, in contrast to these simulations of Paper I, Model 
WOO does not explode. At t » 390 ms the bipolar oscillations 



reach their maximum amplitude. In the further evolution they be- 
come weaker and on average the shock radius decreases (Pig.[5j. 
The slow decay of activity is interrupted by several short phases 
of stronger shock expansion and bipolar oscillation, which occur 
quasi-periodically every 50-100 ms. When we stop the simula- 
tion at t — Is the shock has retreated to a radius of only 70 km 
on average. 

Models W00S, W05S, and W00V, in which a slowly con- 
tracting neutron star was assumed, evolve qualitatively very 
similar to Model WOO. However, with increasing contraction 
timescale the oscillation period becomes longer (up to 100 ms) 
and the growth rate of the low-mode instability decreases. All 
these models are dominated by an I — 1 SASI mode and none of 
them is able to explode. 

Also Model WOOF with its rapidly contracting inner bound- 
ary evolves initially quite similar to Model WOO (Fig. |5). 
However, all timescales are shorter: The oscillation amplitude 
starts to grow already after 50 ms, the shock becomes non- 
spherical at < * 130 ms and convection sets in at t as 160 ms. 
Furthermore, Fig.|6]shows that in this model the / = 2 mode (i.e. 
oscillation between prolate and oblate states) is initially more 
strongly excited than the I = 1 mode, which starts to dominate 
only just before the onset of the explosion. 

In contrast to the models with slower boundary contrac- 
tion, the continuous neutrino heating in Model WOOF is strong 
enough to trigger an explosion at f exp = 194 ms. This difference 
is caused by the fact that the faster contraction leads to gravita- 
tional energy release (the accreted matter heats up by compres- 
sion) and thus to higher neutrino luminosities (see Sec. 16.21 for 
further discussion). The anisotropic gas distribution caused by 
the low-mode oscillations becomes frozen in when the shock 
accelerates outward. The shock develops a prolate deforma- 
tion and a single accretion funnel forms in the northern hemi- 
sphere. Since the explosion attains a large-scale asymmetry, the 
anisotropic distribution of the ejecta exerts a strong gravitational 
force that causes an acceleration of the newly formed neutron 
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Fig. 8. Entropy distribution of Models W12F-C (left column) and W12F (right column) for several times. Model W12F-C quickly 
develops anisotropics because of the onset of convection, whereas in Model W12F convection is initially suppressed and low-mode 
SASI oscillations become visible after about 100 ms. After these oscillations have grown to large amplitude and have begun to 
trigger convection also in Model W12F, the two models explode in a qualitatively very similar way, although the detailed structure 
and asymmetry of the postshock flow and supernova shock are clearly different. (Color figures are available in the online version 
of our paper. ) 



star (see Paper I for details about this process and the proce- 
dure of evaluating (postprocessing) our simulations for the re- 
sulting kick velocity of the neutron stai0) Due to the miniscule 



3 Due to the fact that the neutron star core is replaced in our simu- 
lations by an inner grid boundary and thus anchored at the grid center, 
our conclusions on the neutron star kick velocities need to be confirmed 
by ind ependent hydrodynamic m odels without such a numerical con- 
straint. [BurrgwjTeTal] (2006, 2007) seem to be unable to reproduce our 
findings with their recent simulations, using a code setup that allows 
the neutron star to move and employing a strictly momentum conserv- 
ing implementation of the gravitational effects in the fluid equation of 
motion. However, the numerical consequences of their new treatment 
of the gravity source term, and in particular the supposed superiority 
compared to other (standard) treatments, mu s t be demons trated by de- 
tailed numerical tests, which iBurrows etsU ([2006. 2007) have so far 
not presented. The exact reasons for the potentially discrepant results 
are therefore unclear to us and may be manifold. We strongly e mpha- 
size here that the neutron star kicks reported by Scheck et al. (2006) 
were calculated by two independent methods of postprocessing analy- 
sis. First, making use of total linear momentum conservation, the neu- 
tron star recoil was estimated from the negative value of the momentum 
of all gas on the computational grid at the end of the simulations (note 
that using an inner grid boundary at some radius r > leads to momen- 
tum transfer to the excised inner core so that the gas in the simulation 
domain does not retain zero z-momentum). Second, the different time- 
dependent forces that can contribute to the neutron star acceleration, 



boundary luminosity the energy of the explosion remains rather 
low (0.37 x 10 51 erg at 750 ms after bounce, see Table [1] and 
0.5 x 10 51 erg for the extrapolated value at 1 s), but the neutron 
star attains a fairly high kick velocity (v ns « 200 km/s at 750 ms 
post bounce and estimated 350 km/s for r = Is). 

4.3. Models with typical explosion energies 

While the simulations discussed so far demonstrate clearly the 
existence of a non-radial instability that is not convection, they 
were based on the assumption that the core neutrino luminosities 
are negligibly small. In contrast, in W12F and W12F-C bound- 
ary luminosities were assumed such that the explosion energies 
reached values close to those considered to be typical of core- 
collapse supernovae. An overview of the evolution these models 
can be obtained from Figures [8] and [9] where we show entropy 
distributions at several times and the mass-shell plots, respec- 
tively. 



i.e., the gravitational force excerted by the anisotropically distributed 
matter around the neutron star as well as the momentum transfer associ- 
ated with anisotropically accreted or outflowing gas of the neutron star, 
were added up and then integrated in time. Both of these completely 
independent approaches led to estimates for the neutron star kick veloc- 
ities in very nice agreement with each other. 
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Fig. 9. Upper panel: Same as Fig. [5] but for Model W12F. 
After an initial phase, in which the model remains nearly spher- 
ically symmetric, the SASI becomes strong enough to deform 
the shock and to trigger convection. This model explodes at 
t 160 ms (marked by the vertical dotted line), after the oxygen- 
enriched silicon shell is has fallen through the shock. Lower 
panel: Mass-shell plot for Model W12F-C. Convective activ- 
ity starts to deform the shock in this model at t « 60 ms, and 
the explosion occurs already at t « 120 ms, before the oxygen- 
enriched silicon layer (whose inner boundary is indicated by the 
line labelled with "O in Si") has fallen through the shock. ( Color 
figures are available in the online version of our paper. ) 

In Model W12F-C, in which large initial seed perturbations 
were assumed (cf. Table [TJ, the first convective bubbles form at 
t » 60 ms, and at t ss 90 ms the whole gain layer has become 
convective (see Fig. [8j. From this time on the total energy in the 
gain layer rises continuously and already at f exp « 120 ms the 
first zones acquire positive total energy and the model explodes. 
The initially weakly perturbed Model W12F behaves differently 
in the first 200 ms. There is no sign of convection and for the 
first 100 ms the shock radius evolves as in a corresponding one- 
dimensional model. However, as in Model WOO a weak / = 1 
oscillation mode is present in the postshock flow already at early 
times (f » 30 ms) and grows exponentially to large amplitudes. 
At about t a \ as 150 ms steep convectively unstable entropy gra- 
dients are generated behind the oscillating shock and within two 
cycle periods a situation develops that strongly resembles model 
W12F-C at the onset of the explosion. Also Model W12F ex- 
plodes, though a bit later than Model W12F-C, at t — 164 ms. 



Fig. 10. Evolution of the explosion energy (thick) and the neu- 
tron star velocity (thin) for Models W12F (solid) and W12F-C 
(dotted). 

Although the pre-explosion evolution and the explosion 
timescales of the two models are different, the models behave 
quite similar after the explosion has set in. The convective struc- 
tures merge and downflows form at the interface between ex- 
panding, neutrino-heated gas and the matter with lower entropy 
just behind the shock. The number of downflows decreases with 
time and from t as 200 ms on a single downdraft dominates the 
anisotropic gas distribution. Its position differs in the two mod- 
els, as does the shape of the shock, but the explosion energies 
and even the neutron star velocities grow nearly in the same way 
after t ~ 0.3 s and reach essentially the same values at the end of 
our simulations (Fig.fTOii. 

5. The linear phase of the instability: identification 
of the AAC 

We now turn to a detailed investigation of the question which 
physical mechanism is responsible for the SASI that we have 
seen in the models discussed in the previous section. 



5.1. Motivation and method 



In a number of studies (e.g., Galletti & Foriizzol 


2005; 


Burrows et al. 


2006; Ohnishi et all l2006t IScheck et al.l 


2006; 


Foglizzo et al. 


2007; Yamasaki & Yamada 2007) the advective- 



acoustic cycle was identified or invoked as the cause of the 
SASI oscillati ons that were found in these studies to occur as 
observed by IBlondin et al.l (l2003 | ). This inte r pretati on is cur- 
rently challenged by iBlondin & Mezzacapp a (2006), who ad- 
vocate as an explanation of the SASI modes a purely acous- 
tic process, which is driven by sound waves traveling solely in 
non-radial direction (Blondi n & Shawl I2007I) . One difficulty of 
deciding about the correct interpretation is due to the fact that 
the oscillation timescale of the SASI can either be understood 
as the acoustic timescale along a well chosen transverse path, 
or the advection time down to a suitably chosen coupling ra- 
dius. From the physics point of view, however, the foundations 
of the advective-acoustic mechanism are well documented (see 
the papers cited in Sect. |2.2| and the references therein), whereas 
the purely acou stic mechanis m is still incompletely under stood 
dLamindl2007l) . In particular [Blondin & Mezzacappa! d2006l) ar- 
gued that the existence of a different gradient of the momentum 
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Fig. 11. Time evolution of the amplitude of the dominant spher- 
ical harmonics mode of the pressure, normalized by the ampli- 
tude of the I = mode, as function of radius for Models WOO, 
WOOF and W12F. The solid lines are the minimum, average, and 
maximum shock radius, the dotted line is the gain radius, the 
dashed line is the neutron star surface (defined as the location 
where the density is 10 1 1 g cm -3 ), and the dash-dotted line marks 
the position, Rv(t), of the largest velocity gradient. A low-mode 
oscillation develops in the postshock flow. A pronounced phase 
shift is visible at a radius R^t) that agrees well with the position 
of the largest velocity gradient. The "noise" (short-wavelength 
sound waves) visible in the early phase after bounce is caused 
by the shock propagation and is not related to the advective- 
acoustic cycle. ( Color figures are available in the online version 
of our paper.) 



flux on both sides of the shock is responsible for the instability. 
This argument, however, is so inconclusive that it was used by 
iNobuta & Hanawal dl994 Fig. 10) in order to reach the exactly 
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Fig. 12. Time evolution of the amplitude of the dominant spher- 
ical harmonics mode of the quantity A(r, t, 6) of Eq. ( fT9] >, dis- 
played as function of radius for Models WOO, WOOF and W12F. 
The lines have the same meaning as in Fig. QT| As in the latter 
figure, a zebra-like pattern becomes visible here already several 
ten milliseconds after bounce a zebra-like, indicating that matter 
with a nonvanishing lateral velocity component is advected from 
the shock towards the neutron star. (Color figures are available 
in the online version of our paper.) 



opposite conclusion, namely the stability of a stationary shock 

in an accretion disk. 

Ind ependent of any timescale consideration, Fog lizzo et al.l 
d2007h were able to directly measure the efficiencies of both 
advective-acoustic and purely acoustic cycles using a WKB ap- 
proximation, i.e. for perturbations whose wavelength is shorter 
than the size of the flow gradients near the shock. For every un- 
stable eigenmode for which this quantitative estimate was pos- 
sible, it showed the stability of the purely acoustic cycle and the 
instability of the advective-acoustic one. The WKB approxima- 
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should correspond to abrupt deceleration while smoothly decel- 
erated flows should be stable. Moreover, the amplification factor 
Q during one oscillation is compared to the value measur ed in 
the simpler setups st udied by Blon din & Mez zacappa ( 2006}) and 
IFoglizzo et all d2007l) . 

5.2. Extracting eigenfrequencies from the simulations 

In Fig. Q~2] advected perturbations are displayed by the ampli- 
tudes of the largest modes of the spherical harmonics of a quan- 
tity A(r, f, ff), which turns out to be particularly useful for a quan- 
titative analysis of the SASI. It is defined as 



1 d 

A(t, r,6) = — — (v e (t, r, 6) sin 6) , 
sm8 88 



(19) 



Fig. 13. Oscillatory growth of the amplitudes ci\ and 02 of the 
I — 1 , 2 spherical harmonics components of the quantity A(r, t, 9) 
of Eq. ( TT9b at the gain radius of Model WOOF. (A color figure is 
available in the online version of our paper. ) 



with r A being the divergence of the lateral velocity compo- 
nent, i.e., A = rdiv(veCfl), which scales with the size of the lat- 
eral velocity of the fluid motion. At the gain radius, its expansion 
in spherical harmonics Yi „,{6, 0) is written as 



tion is unfortunately unable to treat accurately the lowest fre- 
quency modes, whose wavelength is comparable to the radius of 
the shock. This argument in principle leaves room for alternative 
explanations of the instability of the lowest frequency modes. 
This is why we do not discard the possibility of a purely acous- 
tic, unstable cycle a priori, despite its unsatisfactory theoretical 
foundation. 

The quantities and results shown in Figs. [TTl - [T6l in the 
present paper are supposed to characterize the development of 
the SASI in a time-dependent environment and to serve com- 
parison of the SASI properties with the expectations of either 
an advective-acoustic or a purely acoustic process. Using a pro- 
jection of perturbations on spherical harmonics, the time evo- 
lution of the radial structure of the most unstable eigenmode is 
visualized, and the oscillation frequency a> r and growth rate <y, 
can be measured. The oscillation timescale is then compared to 
some reference timescales associated with advection and acous- 
tic waves. The acoustic timescales chosen for this comparison 
are T^ nd , computed along a radial path crossing the shock di- 
ameter and back, and T^' und , computed along the circumference 
at the shock radius (i.e., immediately behind the shock position): 



t™ 1 = 2 

sound 



r Rsh dr r R ' b dr 4R ib 

+2 + — i, (15) 

jR ib c-v J R .. c + v c s , ib 



sound 



JR, 

2jrR sh 



(16) 



Following iFoglizzo et al.l d2007l) . the reference timescales cho- 
sen for the advective-acoustic cycle are the advection time rj dv 
from the shock to the radius Ry of maximum deceleration, and 
an estimate of the full cycle timescale rj ac based on a radial ap- 
proximation for simplicity: 



C R ' h dr f* sh d 
Jfiv M J«v c ~ 



(17) 
(18) 



The consistency of the advective-acoustic interpretation is fur- 
ther tested by comparing the timescale of deceleration |dv/dr|~' 
with the oscillation time of the instability. If velocity gradients 
are indeed responsible for the acoustic feedback, unstable flows 



A(f,/? g (f),0) = J] fl/(f)F/ - o(0 ' O) ' 



(20) 



1=0 



where due to the assumption of axisymmetry only m — has to 
be considered. 

For / > 0, the spherical harmonics coefficients a/ of 
this quantity a re proportional to t he ones of the shock dis- 
placement (see IFoglizzo et al]|2006l Appendix F), so A(t,R g ,9) 
contains basically the same information as R s (t,ff). As 
Blondin & Mezzacappa (2006), we prefer to consider a local 
quantity A(t, R g (t), ff) in the postshock layer here rather than the 
shock dis placement §R = Rs (0) - (R s )e (used in iBlondin et al.l 
2003 and Ohni shi et al.l l2006). because A is much less affected 
by noise (A(t) = for a non-stationary spherical flow, whereas 
R s (t) is varying) and allows one to measure the oscillation pe- 
riod and the growth rate much more sensitively than it is pos- 
sible by using R s . Tests showed that for our models, in which 
relatively large seed perturbations were imposed on the infalling 
stellar matter ahead of the shock, A as defined in Eq. (T% 
yields a cleaner measure of the SASI even for very low am- 
plitudes than the perturbed e ntropy or pressure considered by 
Blondin & Mezzacappa (2006). As an example, the absolute val- 
ues of the coefficients a\ and ai are shown as functions of time 
for Model WOOF in Fig. Qj] 

For a given mode / the oscillation period r 0SCj / can be deter- 
mined from the minima of \ai(t)\, which occur at times t — t " . , 

1 ' v / " mm,/ 

when « is a counter for the minima. The detection of the minima 
works reliably only when the amplitude is large enough (it there- 
fore fails in the first 10-20 ms) and is also not feasible when con- 
vective instabilities involve a broad range of frequencies in the 
nonlinear phase. During one cycle of mode I the corresponding 
coefficient a/(f) becomes zero twice, therefore the cycle period 
can be measured as 



T osc,z(Cin,z) _ f I 



77+1 _ 

min,/ min,/ 



(21) 



The evolution of the period of the I = 1 modes, r 0SCj i, evaluated 
with Eq. (|2TT >. is displayed for three of our models in Fig. [15] 

In order to measure the cycle efficiency, Q, we use again the 
coefficients a/ defined in Eq. (f2Qb . We detect the positions of the 
maxima of \ai(f)\, which occur at times t = f x m . dxl (« now being 
the counter for the maxima): if the oscillations of mode / are 
dominated by the (k + 1 )-th harmonic, \ai(f)\ has 2(k + 1 ) maxima 
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during one fundamental cycle period T l aac , so the amplification 
per fundamental cycle can be measured as 



!2(Cx./) = exp 



OJj 





fl,(f" +1 ,)" 

,v max,F 




aiif- 1 ,) 

. ' v max,/ 7 J 



k+\ 



(22) 



This method fails, if several of the harmonics are excited with 
similar strength. However, these phases can be identified and 
typically one of the harmonics dominates clearly (mostly the 
fundamental mode, k = 0). The efficiencies measured by using 
Eq. d22"l > are shown in Fig. [171 

5.3. Interpretation of the results 

The projection of acoustic and advected perturbations on spher- 
ical harmonics reveals that the shock oscillations are associated 
with coherent pressure fluctuations and with the downward ad- 
vection of perturbations produced at the shock. This associa- 
tion is visible in all simulated cases, and is clearly illustrated 
by Figs.[II]and[T2]for Models WOO, WOOF, and W12F. 

The pattern of the pressure perturbations in Fig. [TT] reveals 
the presence of a particular radius where a phase shift occurs. 
The dash-dotted line in these figures is defined as the radius Ry 
where the velocity gradient of the unperturbed flow has a local 
extremum. This particular radius seems to have an important in- 
fluence on the properties and behavior of pressure perturbations; 
in all studied cases the two radii coincide: 



(23) 



This striking coincidence might be interpreted as the conse- 
quence of a particularly efficient coupling between advected 
and acoustic perturbations in layers where the accretion flow is 
strongly decelerated. In order to test this hypothesis, we have 
compared the wavelength 2nv/o> r (to r was defined in the context 
of Eq. [Tol l of advected perturbations at radius Ry to the length 
scale |d In v/drp 1 of this deceleration zone. An efficient coupling 
is expected if the flow velocity varies on scales shorter than the 
wavelength of advected perturbations: 



din v 



dr 



2nv 



(24) 



This hypothesis is confirmed by our simulations for Models 
W12F and WOOF in Fig. [14] In the case of the neutrinoless 
Model W00FA, where only a very weak SASI mode develops, a 
prominent deceleration peak is indeed absent. Therefore the in- 
terpretation of /?v within the framework of the AAC as an effec- 
tive coupling radius between advected perturbations and acous- 
tic feedback, i.e. 



R c = Ry ~ Ra 



(25) 



is consistent with the results of our simulations: efficient cou- 
pling of advective and acoustic perturbations requires a well lo- 
calized deceleration peak, in which case a strong SASI mode can 
be expected to develop. We would like to mention that the ex- 
trema of the flow deceleration that are present in our simulations 
with approximative neutrino transport are also found in simula- 
tions with more sophisticated energy-dependent neutrino trans- 
port and therefore seem to be generic features of the neutrino- 
cooled accretion flow. 

It is interesting to compare the oscillation period t osc mea- 
sured for our models with the timescale rj^ of the fundamental 
AAC mode, approximated by the advection time t J dv of the fluid 
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Fig. 14. Absolute values of the radial derivative of the radial 
velocity component as functions of radius for Models W12F, 
WOOF, and W00FA at several times. The gray-shaded area indi- 
cates the range of values of t^ c during these times. In the mod- 
els including neutrinos a pronounced "deceleration peak" forms 
with a maximum value significantly higher than t~^ c , whereas 
such a feature is absent in Model W00FA. 



and Fig. |T3T >. A systematic comparison between the measured 
oscillation timescale, the advection timescale, and the acoustic 

j-ad „„a -Jat 



timescales r . and 

sound 



T sound ^ S snown m Fig- CLS f° r s i x °f our 



moving from the shock to Ry (r^ 



adv' 



eight models. In all models except W12F, the advection time 
is very close to the oscillation period, whereas in Model W12F 
we find T^ d < t 05C . In the light of the perturbative analysis of 
Eqs. [T7] [H |Yamasakf& Yamada (2007), this finding can be interpreted as 
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Fig. 15. Evolution of the I = 1 mode oscillation period, r osc , the 
advection time rj dv from the shock to the radius /?v of strongest 
deceleration, and the time rj ac for a radial advective-acoustic cy- 
cle, for Models WOO, W12F and W05 V. The advection time rj dv 
agrees with the oscillation time for Models WOO and W05V. In 
Model W12F, the oscillation period is longer than both tL and 
Ta ac , which can be explained by the consequences of strong neu- 
trino heating (see text). The oscillation period is also compared 
to the radial sound travel time r™L d through the shock cavity and 

back, and the maximum lateral sound travel time r lat , just be- 
sound j 

hind the shock. r osc is close to d for Model WOO, but is closer 
to r lat , for Model W05V. The vertical dotted lines enclose the 

sound 

time intervals considered for the evaluations of Figs.[16land[17l 
a consequence of the strong neutrino heating in Model W12F3. 
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Fig. 16. Advection time tL of the fluid moving from the shock 
to the radius of maximum deceleration, Ry, and the two acoustic 
times T™' und (lower grey symbols) and r lat . (upper grey sym- 
bols), versus the oscillation period, t osc , for Models WOOF, WOO, 
WOOS, W05S, W05V and W12F. The data used in this figure 
are selected from phases in which the oscillations can be clearly 
identified and in which the flow is quasi-stationary (see Fig. IT3Tl. 
While r lat , is clearly too long in all cases and r lad , is too 

sound J ° sound 

short for all models except WOO and WOOF, the oscillation pe- 
riod is well approximated by tJ. for all models except W12F. 
The special role of Model W12F can be explained by the effects 
of strong neutrino heating (see text). 



Yamasaki & Y amadal (120071) measured the continuous transition 
of the eigenfrequency from the oscillatory SASI to the purely 
growing (a> r = 0) convective instability when neutrino heating 
is increased. According to their work, the oscillation frequency 
a> r of the SASI is sensitively decreased by the effect of buoy- 
ancy in the gain region, result ing in a significantly longer oscil- 
lation timescale (see Fig. 3 in Yamasa ki & Yama da 2007). This 
agrees well with our results, comparing in particular Models 
WOOF and W12F, whose prescribed contraction of the lower ra- 
dial grid boundary is similar, but the latter model has a much 
larger core (and higher total) neutrino luminosity (see Table [T), 
much stronger neutrino heating, stronger buoyancy, and there- 
fore a larger value of t osc - In contrast, the advection timescale 
is increased by convection to a lesser extent (see Fig. 4 in 
Yamasaki & Yamada 2007), consistent with our finding of the 
data points (t osc ,tL ) for Model W12F lying below the diago- 
nal, dotted line in Fig. [16] The effect of buoyancy can be seen in 
the pressure evolution of Model W12F, shown in the lower plot 
of Fig.QT] where a phase shift <p takes place in the vicinity of the 
gain radius (cf. Eq.fTTIi. 

The effect of buoyancy in Model W12F is also visible in 
Fig. [17] showing the amplification factor Q for six of our eight 
simulated models . The amplification factor has modest values 
between 1 and 3 in most cases, whereas it is as high as Q 6 in 



4 Note that in Model W12F a larger core neutrino luminosity was as- 
sumed than in the other models (see Table[T](. Moreover, the prescribed 



rapid contraction of the nascent neutron star leads to a large accretion 
luminosity. Both contributions to the neutrino emission cause a partic- 
ularly strong neutrino energy deposition in the gain layer. 
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radius of the standing accretion shock in these models is larger 
than in the other cases and therefore the accretion velocities 
in the postshock layer are smaller. This enhances the discrep- 
ancy between the advection time and the radial sound crossing 
time in these models. Given the lack of any better suggestions 
for a unique definition of the timescale of the acoustic cycle 
than the lower and upper bounds considered here, and because 
of the remarkable correlation between the oscillation time and 
T adv' we mter P re t Fig. [16] as a clear support of our hypothesis 
that the SASI oscillations are a consequence of the AAC and 
not of a purely acoustic amp lification process as suggested by 
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Fig. 17. Cycle efficiency, \Q\, as a function of the oscillation pe- 
riod, r osc , for Models WOOF, WOO, WOOS, W05S, W05V and 
W12F. The particularly high values for Model W12F can be ex- 
plained as a consequence of strong neutrino heating in the gain 
layer (see text). 



Model W12F. This high value of Q can be understood as a direct 
consequence of the small value of the oscillation frequency a> r , 
see Eq. ( l22l . because according to Yamasa ki & Yama da (2007) 
stronger neutrino heating sensitively increases r osc , i.e. reduces 
oj r = 27t/t osc , but hardly af fects the growth rate w, (cf. Fig. 2 in 
Yamasaki & Yamada 2007), which appears in the numerator of 
the exponent in Eq. (l22l . Model WOOF exhibits similar trends of 
T adv < T ° sc m ^ 2- ennancemen t as Model W12F, however much 
less strongly. Although in this model the rapid contraction of the 
inner grid boundary leads to a significant accretion luminosity, 
only a very small neutrino flux from the excised inner core was 
assumed and therefore the neutrino heating in the gain layer is 
less strong than in Model W12F. 

We wish to point out that our calculation of the amplification 
factor Q does not rely on any interpretation of the underlying 
mechanism. Interestingly, however, the values between 1 and 
3 are consistent with those measured by Foglizzo et al. (2007, 
Fig. 17) for a shock radius R s h ~ 2Ry in a much simpler con- 
text. From the point of view of the underlying mechanism, these 
values for Q are consistent with numbers obtained by downward 
extrapolation of the efficiency |(2Iwkb of the advective-acoustic 
cycle from the regi on of its validity at larger shock radii (also 
shown in Fig. 17 of lFoglizzo et al.[ 2007). For each of the models 
depicted in Fig. [17] the amplitude of the spread of amplification 
factors can receive a natural explanation in the context of the 
advective-acoustic mechansim: the contribution of the acoustic 
cycle can be either constructive or destructive, depending on the 
relative phase of the two cycles, which varies with time as the 
size of the cavity evolves. This dispersion can be interpreted as 
a measure of the efficiency of the acoustic cycle, which is con- 
sistently smaller than unity. 

The comparison of the oscillation periods with the acous- 
tic timescales shows that t osc is similar to the radial acoustic 
timescale r™^ ind only in Models WOOF and WOO. It is longer by 
up to 30% in the case of Model WOOS and by up to a factor 
of about two in the case of Models W05S and W05 V (Figs. [15] 
and [T6ll. For all models, the upper bound of the acoustic time, 
T sound' * s a l wavs larger than t osc by 20-50%. Note that the setup 
of Models W05S and W05V (with slow contraction of the in- 
ner grid boundary and non-negligible core neutrino luminosities 
and thus significant neutrino heating) was chosen such that the 



The flow properties that are consistent with an advective- 
acoustic cycle as the physical mechanism for the SASI are sum- 
marized as follows: 

(i) The acoustic structure of the unstable modes is strongly cor- 
related with the structure of the velocity gradients (Fig.[TTI). 

(ii) the deceleration region is more localized in unstable flows, 
while smoothly decelerated flows are more stable (Fig[T4li. 

(iii) the advection time rj dv is in good agreement with the oscil- 
lation period of the instability (Fig. [Toll. 

(iv) the typical efficiencies Q ~ 1-3 computed in Fig. [T7] are 
consistent with those extrapolated from the WK B analysis 
of the advective-acoustic cycle in Fig. 17 of lFoglizzo et al.l 
(2007). Their dispersion smaller than unity is consistent with 
the expected marginal effect of the purely acoustic cycle. 

In contrast, a purely acoustic interpretation would have to 
consider the properties (i), (ii), and (iii) as coincidences, and the 
distribution of efficiencies Q remain uninterpreted. Although the 
oscillation period is consistent with an acoustic timescale along a 
carefully chosen acoustic path, Fig. [16] indicates that this acous- 
tic path should be close to radial for some models, and much 
more lateral in others. A theory of the purely acoustic instability 
would have to explain this behavior. 

Without claiming that our present knowledge of the 
advective-acoustic theory is fully satisfactory in the complex 
core-collapse context, its mechanism is understandable from the 
physics point of view and allows us to explain several features 
of the simulations, which would not be understood otherwise. 



6. Interpretation of the nonlinear phase 

In the following we will discuss our simulations during the non- 
linear phase of the evolution in which the SASI cannot be con- 
sidered as a small perturbation. In particular, we will analyse the 
relation between the SASI and convective instability, as well as 
the role these instabilities play for the explosion mechanism and 
the resulting energy of the explosion. 



6.1. The SASI as trigger of convective overturn 

In models with low core neutrino luminosity convective activ- 
ity does initially not occur because the corresponding instability 
is suppressed in the accretion flow of the neutrino-heating layer 
according to Eqs. ([6]) and ([8). The first large-scale non-radial 
perturbations in the postshock flow of such models are therefore 
caused by SASI oscillations. Once large average lateral veloci- 
ties around 10 9 cms -1 or more are reached in the gain layer at 



18 



Scheck et al.: Multidimensional supernova simulations 



t > fjj (cf. Table Q]), however, also the smaller-scale mushroom- 
like structures that are typical of the onset of Rayleigh-Taylor in- 
stability start to grow. Within only a few more oscillation cycles, 
plumes of neutrino-heated matter and supersonic downdrafts of 
low-entropy matter develop and violent convective overturn sets 
in very similar to what we found in the case of the models de- 
scribed in Paper I. There are two effects that are mainly respon- 
sible for the corresponding change of the flow character; these 
are linked to the unsteady motion and the growing deformation 
of the shock, respectively. 

Firstly, in course of radial expansion and contraction phases 
the shock reaches velocities of <9(10 9 cm/s), which is a signif- 
icant fraction of the preshock velocity. Since the postshock en- 
tropy depends on the preshock velocity in the frame of the shock, 
such fast shock oscillations cause strong variations of the en- 
tropy in the downstream region. Rapid outward motion of the 
shock produces high entropies in the postshock flow, whereas 
phases in which the shock retreats lead to lower postshock en- 
tropies. Periodic shock expansion and contraction thus results in 
alternating layers with high and low entropies, which are com- 
pressed as the accreted matter is advected towards the neutron 
star. With increasing amplitude of the shock oscillations the con- 
vectively unstable entropy gradients between these layers even- 
tually become so steep that the growth timescale of Rayleigh- 
Taylor instabilities shrinks to about 1 ms, which is much shorter 
than the advection timescale. Therefore non-radial perturbations 
are able to grow quickly at the entropy interfaces and vortices 
and mushroom-like structures begin to form (see Fig. [7j>. 
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Fig. 18. Evolution of the sum of the v e and v e luminosities at the 
neutrinosphere and at r — 500 km for Models WOO and WOOF. 
The luminosities decay with time because the largest mass accre- 
tion rates are present at early times (t < 0.2 s). Model WOOF has 
the higher luminosity, because the neutron star contracts faster, 
setting free more gravitational energy. 

Secondly, also the off-center displacement of the accretion 
shock by the I = 1 SASI mode and the shock deformation 
caused by / > 2 modes play an important role when the ampli- 
tudes become large enough. The radial preshock flow hits the de- 
formed or displaced shock at an oblique angle. Since the velocity 
component tangential to the shock is not changed when the gas 
passes through the shock, in contrast to the normal component, 
which is strongly reduced, the flow is deflected and attains a sub- 
stantial lateral velocity component, whose size and sign changes 
during the cycle period, see Fig. [6] As long as the cycle ampli- 
tude is small, the lateral velocity components are also small, and 



the postshock flow remains approximately radial. In the case of a 
strongly deformed shock, however, the postshock flow becomes 
mainly non-radial because the lateral velocity reaches a signif- 
icant fraction of the preshock velocity (up to several 10 9 cm/s, 
i.e. the lateral flow becomes supersonic). For an / = 1 mode, the 
highest negative lateral velocities are obtained when the shock 
has its maximum displacement in the negative z-direction, see 
Fig. |6l upper left panel. A shell of matter with high negative lat- 
eral velocity formed in this phase is advected towards the neu- 
tron star, and half an oscillation period later the highest posi- 
tive lateral velocities are generated right behind the shock when 
the shock expands into positive z-direction (Fig. |6l middle left 
panel). With increasing oscillation amplitude the shock radius 
- and consequently also the advection timescale and the cycle 
period - begin to vary so strongly during one cycle that the 
northern and southern hemispheres run "out of phase" so that 
the shock radii at the north pole (0 = 0°) and at the south pole 
(9 = 180°) reach their maximum values not alternatingly any 
more, but at almost the same time. In this case streams of matter 
with high positive and high negative lateral velocities emerge si- 
multaneously near the north and south pole, respectively. These 
streams collide and one of them is deflected upwards, producing 
a bump bounded by two "kinks" in the shock surface, while the 
other one is directed downwards, forming a supersonic down- 
flow (see Fig. [6] lower left panel and Fig.[8j middle right panel), 
a phenomenon that we have also observ ed in the simu l ations 
of Pa per I and that was also reported by iBurrows et al.l (2006, 
l2007h . 

Large-amplitude SASI oscillations are thus able to trig- 
ger nonlinear convective overturn even in models in which the 
growth of buoyancy instabilities is initially suppressed because 
of unfavorable conditions in the accretion flow as discussed in 
the context of Eqs. (|6]l and ((SJ. 



6.2. From SASI oscillations to explosions 

Why is Model WOOF able to develop an explosion while Model 
WOO and the models with even slower boundary contraction 
(W00S, W05S, W05V) do not explode? Models WOO and WOOF 
differ in the assumed contraction of the nascent neutron star, 
i.e. in the parameters describing their final inner boundary ra- 
dius, R f ib , and the contraction timescale, to,. A smaller value of 
R f ib implies that the matter accreted on the forming neutron star 
sinks deeper into the gravitational potential and thus more grav- 
itational energy is released. The smaller value of f;b causes this 
release of energy to happen earlier. Most of the potential en- 
ergy that is converted to internal energy by pd V-work is radiated 
away in the form of neutrinos. Consequently, the neutrino lumi- 
nosity that leads to heating in the gain layer is much higher at 
early times in the case of Model WOOF (Fig. [181. 

Yet, these high luminosities alone are not sufficient to start an 
explosion. This is demonstrated by a one-dimensional simula- 
tion with the same boundary parameters as Model WOOF, which 
fails to explode. It is well known that in the multi-dimensional 
case convection leads to an enhancement of the efficiency of 
neutrino energy deposition in the gain layer, on the one hand 
because non-radial convective motions stretch the time fluid el- 
ements can stay in the gain layer and are thus exposed to ef- 
ficient neutrino heating in the vicinity of the gain radius, on 
the other hand because high-entropy, neutrino-heated matter be- 
comes buoyant, expands quickly, and thus cools, which reduces 
the energy loss by the reemission of neutrinos. The former of 
these two effects effectively leads to an increase of the advec- 
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Fig. 19. Evolution of mass in the gain layer, 
AMgain, advection timescale T a d v from the 
shock to the gain radius, and neutrino heat- 
ing rate integrated over the whole gain layer, 
Q gmn , for Models WOO and WOOF (solid 
lines). The dotted lines are results from cor- 
responding ID simulations. When the SASI 
becomes nonlinear (f » 350 ms for WOO, 
t « 150 ms for WOOF) and triggers convec- 
tion (~30 ms later), T adv and AM gilin begin to 
grow. As a consequence, Qgafa increases sig- 
nificantly compared to the ID simulations. In 
the case of Model WOOF the enhanced heat- 
ing is strong enough to gravitationally unbind 
the matter in the gain layer and to drive an 
explosion. The advection timescale shows an 
increase also in the ID counterpart of Model 
WOOF at / a 150 ms, because a composition 
interface of the progenitor star falls through 
the shock at this time and the strong decrease 
of the mass accretion rate leads to a tran- 
sient expansion of the shock and of the gain 
layer. Nevertheless the ID model does not de- 
velop an explosion because without the aid 
of multi-dimensional effects neutrino heating 
cannot become powerful enough. In Model 
WOO the neutrino energy deposition rate is so 
low that even an increase by almost a factor 
of two between 350 ms and 400 ms is not suf- 
ficient for an explosion. 



tion timescal e of accreted matte r from the shock to the gain ra- 
dius (see also lBuras et al.l2006bl) . as a consequence of which the 
mass in the gain layer becomes larger. The same effect can also 
be produced by large-amplitude SASI oscillations, because such 
non-radial motions expand the average shock radius, thus lead- 
ing to smaller postshock velocities, and deflect the postshock 
flow in non-radial direction, also leading to a longer advection 
time of accreted matter through the gain layer. 

In Model WOOF we observe such a rise of the advection 
timescale starting at t » 150 ms (Fig. [T9l> when the postshock 
flow becomes strongly non-radial, but violent convective over- 
turn has not yet set in (Fig. [6] right middle panel). This increase 
of the advection time leads to a significant growth of the inte- 
grated neutrino heating rate in the gain layer, an effect that be- 
comes even more pronounced when the convective activity gains 
strength (f > 170 ms). Initially the total specific energy of most 
of the matter in the gain region is in a narrow range around 
- 1 1 Me V per nucleon, but the distribution of specific particle en- 
ergies becomes broader by the influence of the large-amplitude 
SASI and of convective overturn (Fig. 120b. Due to the large en- 
ergy deposition by neutrinos the mean value of the total energy 
rises and ultimately some fraction of the matter in the gain layer 
acquires positive total energy and the explosion sets in. Also in 
Model WOO we see enhanced neutrino heating (up to two times 
higher than in the corresponding one-dimensional simulation) 
from t * 350 ms on, caused by a combination of nonlinear SASI 
motions and convective activity (Fig. [T9| >. However, due to the 
low accretion rate at this late time the neutrino luminosity and 
thus the neutrino heating rate are much lower than in Model 
WOOF at t » 200 ms. The total energy in the gain layer of Model 
WOO increases only temporarily by about 1 MeV per nucleon but 
then drops again soon and continues to decrease slowly later on 
(Fig. [20b . The distribution of specific energies of matter in the 



gain layer does not become very broad and none of the matter 
gets unbound. In both the Models WOO and WOOF the specific 
kinetic energy in the gain layer remains relatively small (only 
about 1 MeV/nucleon, see Fig.l20l>. 

Different from iBlondin et all (120031) we do not observe a 
continuous increase of the kinetic energy associated with lateral 
(turbulent) motion of the matter b ehind the shock. In their sim- 
ulations without neutrino effects. IBlondin et"aLl f2003) observed 
that the SASI oscillation can redistribute some of the gravita- 
tional binding energy of the accreted matter from coherent fluid 
motion to turbulent energy, in fact with sufficient efficiency to 
drive an expansion of the accretion shock. Since some of the tur- 
bulent material h ad obtained positive total energy at the end of 
their simulations, Blond irTet al.l d2003) concluded that the SASI 
in their calculations was able to lead to an explosion. We do not 
see this kind of process going on in our simulatio ns (in agree- 
ment with the results of Burr ows et al.ll2006l 120071) . The reason 
for this discrepancy may be the inclusion of neutrino heating 
and cooling in our models. It is possible that the energy loss by 
neutrino emission below the gain radius prevents the efficient 
conversion of gravitational binding energy to turbulent energy. 
Another reason may be the use of different c onditions at the outer 
radial grid boundary in our models; while iBlondin et all ([2003 ) 
assumed steady-state accretion and thus held the mass accretion 
rate fixed with time, the stellar progenitor structure employed in 
our work leads to a continuous decrease of the mass accretion 
rate at the shock. Therefore less total kinetic energy is available 
that can be converted to turbulent motions by the distorted ac- 
cretion shock. 

In our simulations a growth of the turbulent kinetic energy 
of the matter in the gain layer is definitely not the reason for 
starting the explosions. The corresponding lateral kinetic en- 
ergy never exceeds a few 10 49 erg in any of our models. This 
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energy per baryon (thick solid line) and kinetic energy per 
baryon (dashed line) versus time in the gain layer of Models 
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interval that contains 90% of the mass of the gain layer. 



is well below the size of neutrino energy deposition and of the 
energy needed for unbinding matter and triggering an explosion. 
Nevertheless, the non-radial flow associated with the SASI is 
certainly helpful, in combination with convection actually cru- 
cial for making the neutrino-heating mechanism work. The fail- 
ure of one-dimensional simulations with the same treatment of 
the neutrino physics clearly demonstrates the importance of non- 
radial fluid instabilities, convection and the SASI, for a success 
of the neutrino-driven explosion mechanism. These hydrody- 
namic instabilities affect the gas motion in the gain layer such 
that the advection timescale from the shock to the gain radius 
is effectively increased. This enhances the efficiency of neutrino 
energy deposition by allowing more matter to be exposed to the 
intense neutrino flux near the gain radius for a longer time. Thus 
both convection and the SASI can be considered as "catalysts" 
that facilitate neutrino-driven explosions rather than being direct 
drivers or energy sources of the explosion. As a consequence, 
explosions in multi-dimensional simulations, i.e. even with the 
support by convective overturn and the SASI, still require the 
presence of strong enough neutrino heating. Our set of simula- 
tions clearly demonstrates that only in the case of a sufficiently 
large neutrino luminosity and thus only for sufficiently powerful 
neutrino heating behind the shock, the models are able to de- 
ve lop an explosion. Th is finding is in agreement with the results 
of Ohnish iet alj d2006l) . who also obtained an explosion by neu- 



trino heating only in a simulation with high neutrino luminosity, 
while lower-luminosity cases failed to explode. 

6.3. The importance of the seed perturbations 

The comparison of our results shows that explosions during the 
first second after core bounce do not only require the neutrino 
luminosities to be large enough but also that the SASI and con- 
vection are able to reach the nonlinear phase sufficiently quickly. 
Whether this is the case or not depends on their growth rates, 
which in turn depend on the properties of the postshock flow. 
The latter are a complex function of the progenitor structure, 
the neutron star contraction, and neutrino cooling and heating 
in the layer between neutron star and shock. Last but not least, 
also the size of the seed perturbations, i.e. the inhomogeneities 
present in the collapsing star, can play a role for the develop- 
ment and growth of non-radial hydrodynamic instabilities after 
core bounce. 

In our Models W12F and W12F-C (as well as in most of the 
recent simulations with multi-energy group neutrino transport by 
Buras et al. 2006a b) the advection time through the gain layer is 
so short and the growth rate of convective instabilities (Eq.[TJ in 
that region so small that the timescale ratio x of Eq. © remains 
below the critical threshold xo for a linear gro wth of globally 
unstab le modes, i.e., X < Xo ~ 3 ac cording to i Foglizzo et al.l 
(2006), see Eq. (|6]l and also Fig. 3 in iBuras et al . (2006b). This 
means that the fast advection of the flow from the shock to 
the gain radius suppresses the growth of convective modes ac- 
cording to linear stability analysis. However, as explained in 
Sect. |2.U in this case buoyancy can nevertheless drive bubble 
rise and convective instability in a nonlinear way if the initial 
density perturbations 5 in matter falling through the shock are 
large enough, i.e. 5 > 5 c „t exp(-^) according to Eq. (0, with 
6 CI ii being typically of the order of some percent (see Eq.[8]). 

The inhomogeneities in the matter upstream of the shock 
originate from seed perturbations in the progenitor star, 
whose size and amplitude are not well known because three- 
dimensional, long-time stellar evolution simulations for full- 
sphere models until the onset of core collapse have not been 
possible so far (see, e.g.. iBazan & Arnettlll998t iMurphv et ail 
12004 1 Young et alj2005HMeakin & Arnettl2006ll2007allbl) . With 
our assumed initial inhomogeneities in the case of Models W12F 
and W12F-C (see Table [T] and Sect. 13. 21 ), the perturbation am- 
plitudes remain below the critical value 6 CI i t in the former case, 
whereas they become larger than this threshold value in the latter 
case (see Fig.|4|. Therefore, as visible in Fig. [8] the fastest grow- 
ing non-radial instability on large scales is the SASI in Model 
W12F, whereas it is convective overturn in Model W12F-C. Only 
because of the growth of SASI modes does Model W12F also 
develop convective activity in the gain layer, which enhances 
the efficiency of neutrino energy deposition and finally leads 
to an explosion also in this case. The crucial role of these non- 
radial instabilities is demonstrated again by a corresponding one- 
dimensional simulation that does not develop an explosion. In 
Model W12F the SASI is a key feature in the multi-dimensional 
evolution, because the development of convective modes is not 
possible in the first place due to the low initial amplitude of per- 
turbations and the insufficient growth of these perturbations in 
the advection flow from the shock to the gain radius. 

Considering Models W12F and W12F-C, however, no notice- 
able memory of the initial source of the low-mode asymmetries 
is retained during the long-time evolution. Although the early 
postbounce evolution of these two models is clearly different and 
the times of the onset of the explosion differ, the global parame- 
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ters of the explosion become very similar (see Table[TJ. Neither 
the explosion energy nor the neutron star mass and kick velocity 
are strongly affected by the different explosion times, because 
the conditions in the infalling stellar core change only on longer 
timescales and the ejecta energy and neutron star recoil build 
up over a much more extended period of time after the launch of 
the explosion (see Paper I). Since the anisotropies develop chaot- 
ically and in a very irregular way during the nonlinear phase, the 
final ejecta morphology is the result of a stochastic process and 
does not depend in a deterministic and characteristic manner on 
the type of non-radial instability that has grown fastest after core 
bounce. It therefore seems unlikely that observational parame- 
ters of supernova explosions are able to provide evidence of the 
initial trigger of the large-scale anisotropies that develop in the 
early stages of the explosion. Future simulations with a more de- 
tailed treatment of the neutrino transport (instead of our approx- 
imative description) and without the use of the inner boundary 
condition of the present models will have to show whether the 
gravitational-wave and neutrino signals carry identifiable finger- 
prints of this important aspect of supernova dynamics. 

7. Summary and conclusions 

We performed a set of two-dimensional hydrodynamic simu- 
lations with approximative neutrino transport to investigate the 
role of non-convective instabilities in supernova explosions. As 
initial data we used a postbounce model of a 15 M Q progenitor 
star, which had been evolved through core collapse and bounce 
in a computation with detailed, energy-dependent neutrino trans- 
port. For following the subsequent, long-time evolution, the neu- 
tron star core (above a neutrino optical depth of about 100) was 
excised and replaced by a contracting Lagrangian inner bound- 
ary that was intended to mimic the behavior of the shrinking, 
nascent neutron star. The models of our set differed in the choice 
of the neutrino luminosities assumed to be radiated by the ex- 
cised core, in the prescribed speed and final radius of the con- 
traction of the neutron star, and in the initial velocity perturba- 
tions imposed on the ID collapse model after bounce. 

Our hydrodynamic simulations indeed provide evidence — 
supporting previous linear analysis dFoglizzo et al.ll2006l [2007; 
Yamasaki & Yamada 2007) — that tw o different hydrodyn amic 
instabilities, convection and the SASI dBlondin et al.ll2003l) . oc- 
cur at conditions present during the accretion phase of the stalled 
shock in collapsing stellar cores and lead to large-scale, low- 
mode asymmetries. These non-radial instabilities can clearly be 
distinguished in the simulations by their growth behavior, lo- 
cation of development, and spatial structure. While convective 
activity grows in a non-oscillatory way and its onset can be rec- 
ognized from characteristic mushroom-type structures appear- 
ing first in regions with steep negative entropy gradients, the 
SASI starts in an oscillatory manner, encompasses the whole 
postshock layer, and leads to low-mode shock deformation and 
displacement. 

As discussed by Fogli zzo et alJd2006l) . the growth of convec- 
tion is suppressed in the accretion flow because of the rapid in- 
fall of the matter from the shock to the gain radius, unless either 
neutrino heating is so strong and therefore the entropy gradient 
becomes so steep that the advection-to-growth timescale ratio 
(X of Eq. [5j exceeds the critical value xo ~ 3, or, alternatively, 
sufficiently large density perturbations in the accretion flow (cf. 
Eq.[9]) cause buoyant bubbles rising in the infalling matter. While 
convective instability is damped by faster infall of accreted mat- 
ter, the growth rate of the SASI increases when the advection 
timescale is shorter. Moreover, the quasi-periodic shock expan- 



sion and contraction with growing amplitude due to the SASI 
produce strong entropy variations in the postshock flow, which 
can then drive convective instability. Even when the neutrino- 
heated layer is not unstable to convection in the first place, the 
perturbations caused by the SASI oscillations can thus be the 
trigger of "secondary" convection. 

Our detailed analysis of the evolution of the SASI modes in 
our simulations, of their dependence on the model parameters, 
and of the cooperation between convection and the SASI in the 
nonlinear regime revealed the following facts: 

1 . When the SASI reaches large amplitudes and supersonic lat- 
eral velocities occur in the postshock flow, sheets with very 
steep unstable entropy gradients are formed. As a conse- 
quence, low-mode convective overturn grows in a highly 
nonlinear way. Since the SASI and strong convective activ- 
ity push the accretion shock to larger radii, they reduce the 
infall velocity in the postshock layer. Moreover, the flow in 
the neutrino-heating region develops large non-radial veloc- 
ity components and therefore the accreted matter stays in the 
gain layer for a longer time. This increases the energy de- 
position by neutrinos in this region and facilitates the explo- 
sion. However, like convection the SASI does not guaran- 
tee an explosion on the timescale consi dered in our simula- 
tions (in agreement with the findings of lBurrows et alj|2006l 
2007): the kinetic energy associated with the SASI remains 
negligible for the explosion energetics. Therefore sufficiently 
strong neutrino heating and consequently a sufficiently large 
neutrino luminosity are still necessary to obtain explosions. 

2. The growth rate (and amplification) and oscillation fre- 
quency of the SASI depend sensitively on the advection time 
from the shock to the coupling region and the structure of the 
flow in this region. The latter, in turn, depends on the neutron 
star contraction (which has a strong influence on the shock 
radius), on the neutrino luminosities and the corresponding 
heating and cooling, and on the mass accretion rate of the 
stalled shock in the collapsing star. For a wide range of in- 
vestigated conditions (changing the parameters of our mod- 
els), we found the SASI being able to develop large ampli- 
tudes on timescales relevant for the explosion. Therefore the 
SASI turned out to create large-scale anisotropy also in cases 
where convective activity was not able to set in in the first 
place. 

3. In our simulations we could clearly identify a faster contrac- 
tion and a smaller radius of the excised core of the nascent 
neutron star as helpful for an explosion. We therefore con- 
clude that a softer high-density equation of state and general 
relativity, which both lead to a more compact neutron star, 
are favorable for an explosion. This is so because on the one 
hand the accretion luminosities of neutrinos become higher, 
correlated with stronger neutrino heating behind the stalled 
shock, and on the other hand the amplification of the SASI 
increases with enhanced neutrino heating. Since the initial 
growth of convection is damped or suppressed by more rapid 
infall, the presence of the SASI instability and its ability to 
trigger convection as a secondary phenomenon, play a cru- 
cial role for the final success of the delayed neutrino-heating 
mechanism. 

4. The amplitude of the initial seed perturbations in the collaps- 
ing core of the progenitor star, which is not well constrained 
due to the nonexistence of fully consistent and long-evolved 
three-dimensional stellar evolution models, has an influence 
on the question whether convection or the SASI develop 
more rapidly after core bounce. While this can determine 
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how fast the explosion sets in, we found that once the two 
non-radial instabilities are simultaneously present and coop- 
erate in the nonlinear regime, the global properties of the ex- 
plosions are essentially insensitive to the initial phase. Since 
the final anisotropic distribution of the ejecta is the result of a 
very stochastic and chaotic process, it has also lost the mem- 
ory of the early evolution. Neither the explosion energy nor 
the neutron star kick velocity are therefore good indicators 
of the initial seed perturbations that existed in the progenitor 
star and of the type of the fastest growing non-radial insta- 
bility. Future supernova models without the approximations 
used in our simulations will have to show whether the neu- 
trino and gravitational-wave signals, which will be measur- 
able in detail from a galactic event, carry any useful infor- 
mation about this crucial aspect of the postbounce explosion 
dynamics. 

5. While most of the above conclusions and the corresponding 
evaluation of our simulations are independent of the exact 
physical mechanism that is responsible for the growth of the 
SASI, we nevertheless tried to explore this important ques- 
tion, which is still controversially discussed in the literature. 
To this end we analysed our models in the linear regime of 
the SASI and compared the results on the one hand to pre- 
dictions based on the hypothesis that the SASI growth is due 
to an advective-acoustic cycle (AAC), and on the other hand 
to the possibility that the SASI is driven by a purely acoustic 
mechanism. Our analysis shows that many of the observed 
SASI properties are consistent and can be well understood 
with the AAC hypothesis. This is the case for: (1) The os- 
cillation period of the SASI, which agrees well with the ad- 
vection time of perturbations from the shock to a radius /?v 
where the deceleration of the accretion flow is strongest. This 
radius is located in the neutrino-cooling layer somewhat out- 
side of the neutrinosphere. (2) The acoustic structure of the 
unstable modes is strongly correlated with the velocity gra- 
dient in the postshock layer and more SASTunstable flows 
are obtained in more abruptly decelerated accretion flows. 
(3) The amplification factors found for the SASI agree with 
extrapolated WKB results for the AAC. Moreover, the effect 
that stronger neutrino heating causes a larger SASI amplifi- 
cation efficiency can be explained on grounds of an assumed 
AAC. In contrast, our measured oscillation timescales for the 
SASI are not consistent with a uniquely chosen path for the 
sound wave propagation through the shock cavity in all mod- 
els, but would require that the acoustic waves travel more ra- 
dially in some cases and more in angular direction in other 
models. This as well as the other mentioned features are not 
satisfactorily accounted for by the existing theory of a purely 
acoustic instability. We therefore think that our analysis pro- 
vides significant support for the suggestion that the ampli- 
fication of the SASI happens through an AAC rather than a 
purely acoustic mechanism. 

While the presented simulations employed a number of ap- 
proximations like the neutrino transport scheme and the use 
of an inner boundary condition instead of following the evo- 
lution of the neutron star core, we are confident that our main 
findings do not depend on these simplifications. In fact, recent 
long-time stellar core-collapse simulations with detailed, multi- 
energy-group neutrino transport and fully consistent consider- 
ation of the central part of the nascent neutron star basically 
confirm the importance of the SASI and its nonlinear interac- 
tion with convective instability for the viability of the delay ed 
neutrino-driven explosion mechanism (Marek & Jankal l200"7l) . 
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